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Two gapped quantum ground states in the same phase are connected by an adiabatic evolution 
which gives rise to a local unitary transformation that maps between the states. On the other 
hand, gapped ground states remain within the same phase under local unitary transformations. 
Therefore, local unitary transformations define an equivalence relation and the equivalence classes 
are the universality classes that define the different phases for gapped quantum systems. Since 
local unitary transformations can remove local entanglement, the above equivalence/universality 
classes correspond to pattern of long range entanglement, which is the essence of topological order. 
The local unitary transformation also allows us to define a wave function renormalization scheme, 
under which a wave function can flow to a simpler one within the same equivalence/universality 
class. Using such a setup, we find conditions on the possible fixed-point wave functions where the 
local unitary transformations have finite dimensions. The solutions of the conditions allow us to 
classify this type of topological orders, which generalize the string-net classification of topological 
orders. We also describe an algorithm of wave function renormalization induced by local unitary 
transformations. The algorithm allows us to calculate the fiow of tensor-product wave functions 
which are not at the fixed points. This will allow us to calculate topological orders as well as 
symmetry breaking orders in a generic tensor-product state. 
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I. INTRODUCTION - NEW STATES BEYOND 
SYMMETRY BREAKING 

According to the principle of emergence, the rich prop- 
erties and the many different forms of materials originate 
from the different ways in which the atoms are ordered 
in the materials. Landau symmetry-breaking theory pro- 
vides a general understanding of those different orders 
and resulting rich states of matter. ^'^ It points out that 
different orders really correspond to different symmetries 
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in the organizations of the constituent atoms. As a mate- 
rial changes from one order to another order (i.e., as the 
material undergoes a phase transition), what happens 
is that the symmetry of the organization of the atoms 
changes. 

For a long time, we believed that Landau symmetry- 
breaking theory describes all possible orders in materials, 
and all possible (continuous) phase transitions. However, 
in last twenty years, it has become more and more clear 
that Landau symmetry-breaking theory does not describe 
all possible orders. After the discovery of high super- 
conductors in 1986,'^ some theorists believed that quan- 
tum spin liquids play a key role in understanding high 
Tc superconductors'* and started to introduce various 
spin liquids. Despite the success of Landau symmetry- 
breaking theory in describing all kinds of states, the the- 
ory cannot explain and does not even allow the existence 
of spin liquids. This leads many theorists to doubt the 
very existence of spin liquids. In 1987, in an attempt to 
explain high temperature superconductivity, an infrared 
stable spin liquid - chiral spin state was discovered, ^°'*^ 
which was shown to be perturbatively stable and exist as 
quantum phase of matter (at least in a large N limit). 
At first, not believing Landau symmetry-breaking the- 
ory fails to describe spin liquids, people still wanted to 
use symmetry-breaking to describe the chiral spin state. 
They identified the chiral spin state as a state that breaks 
the time reversal and parity symmetries, but not the spin 
rotation symmetry.^* However, it was quickly realized 
that there are many different chiral spin states that have 
exactly the same symmetry, so symmetry alone is not 
enough to characterize different chiral spin states. This 
means that the chiral spin states contain a new kind of 
order that is beyond symmetry description.*^ This new 
kind of order was named*"^ topological order. 

The key to identify (and define) new orders is to iden- 
tify new universal properties that are beyond the local 
order parameters and long-range correlations used in the 
Landau symmetry breaking theory. Indeed, new quan- 
tum numbers, such as ground state degeneracy*^, the 
non-Abelian Berry's phase of degenerate ground states*'^ 
and edge excitations*^, were introduced to character- 
ize (and define) the different topological orders in chi- 
ral spin states. Recently, it was shown that topological 
orders can also be characterized by topological entan- 
glement entropy.*^'*'' More importantly, those quantities 
were shown to be universal {ie robust against any local 
perturbation of the Hamiltonian) for chiral spin states.*'^ 
The existence of those universal properties establishes the 
existence of topological order in chiral spin states. 

Near the end of 1980's, the existence of chiral 
spin states as a theoretical possibility, as well as 
their many amazing properties, such as fractional 
statistics, *°'** spin-charge separation,*"'** chiral gapless 
edge excitations,*^ were established reliably, at least in 
the large iV-limit introduced in Ref. 26. Even non- 
Abelian chiral spin states can be established reliably in 
the large A^-limit.*** However, it took about 10 years to es- 



tablish the existence of a chiral spin state reliably without 
using large TV-limit (based on an exactly soluble model 
on honeycomb lattice).*^ 

Soon after the introduction of chiral spin states, exper- 
iments indicated that high-temperature superconductors 
do not break the time reversal and parity symmetries. So 
chiral spin states do not describe high-temperature super- 
conductors. Thus the theory of topological order became 
a theory with no experimental realization. However, the 
similarity between chiral spin states and fractional quan- 
tum Hall (FQH) states allows one to use the theory of 
topological order to describe different FQH states. Just 
like chiral spin states, different FQH states all have the 
same symmetry and are beyond the Landau symmetry- 
breaking description. Also like chiral spin states, FQH 
states have ground state degeneracies^* that depend on 
the topology of the space. *'^'^" Those ground state de- 
generacies are shown to be robust against any pertur- 
bations. Thus, the different orders in different quantum 
Hall states can be described by topological orders, and 
the topological order does have experimental realizations. 

The topology dependent ground state degeneracy, that 
signal the presence of topological order, is an amazing 
phenomenon. In FQH states, the correlation of any lo- 
cal operators are short ranged. This seems to imply that 
FQH states are "short sighted" and they cannot know 
the topology of space which is a global and long-distance 
property. However, the fact that ground state degeneracy 
does depend on the topology of space implies that FQH 
states are not "short sighted" and they do find a way 
to know the global and long-distance structure of space. 
So, despite the short-range correlations of any local oper- 
ators, the FQH states must contain certain hidden long- 
range correlation. But what is this hidden long-range 
correlation? This will be one of the main topic of this 
paper. 

Since high Tc superconductors do not break the time 
reversal and parity symmetries, nor any other lattice 
symmetries, some people concentrated on finding spin 
liquids that respect all those symmetries and hoping 
one of those symmetric spin liquids hold the key to un- 
derstand high Tf, superconductors. Between 1987 and 
1992, many symmetric spin liquids were introduced and 
studied. ^"^'^^"^^ The excitations in some of constructed 
spin liquids have a finite energy gap, while in others there 
is no energy gap. Those symmetric spin liquids do not 
break any symmetry and, by definition, are beyond Lan- 
dau symmetry-breaking description. 

By construction, topological order only describes the 
organization of particles or spins in a gapped quantum 
state. So the theory of topological order only applies to 
gapped spin liquids. Indeed, we find that the gapped 
spin liquids do contain nontrivial topological orders^'* 
(as signified by their topology dependent and robust 
ground state degeneracies) and are described by topo- 
logical quantum field theory (such as Z2 gauge theory) 
at low energies. One of the simplest topologically or- 
dered spin liquids is the Z2 spin liquid which was first 
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introduced in 1991.^"^'^ The existence of Z2 spin liquid 
as a theoretical possibility, as well as its many amazing 
properties, such as spin-charge separation, ^'^'^'^ fractional 
mutual statistics,^'' topologically protected ground state 
degeneracy,^'* were established reliably (at least in the 
large iV- limit introduced in Ref. 26). Later, Kitaev in- 
troduced the famous toric code model which establishes 
the existence of the Z2 spin liquid reliably without using 
large A^-limit.^'' The topologically protected degeneracy 
of the Z2 spin liquid was used to perform fault-tolerant 
quantum computation. 

The study of high Tc superconductors also leads to 
many gapless spin liquids. The stability and the exis- 
tence of those gapless spin liquid were in more doubts 
than their gapped counterparts. But careful analysis in 
certain large N limit do suggest that stable gapless spin 
liquids can exist. ^^'^^'^^ If we do believe in the existence 
of gapless spin liquids, then the next question is how to 
describe the orders {ie the organizations of spins) in those 
gapless spin liquids. If gapped quantum state can contain 
new type of orders that are beyond Landau's symmetry 
breaking description, it is natural to expect that gapless 
quantum states can also contain new type of orders. But 
how to show the existence of new orders in gapless states? 

Just like topological order, the key to identify new or- 
ders is to identify new universal properties that are be- 
yond Landau symmetry description. Clearly we can no 
longer use the ground state degeneracy to establish the 
existence of new orders in gapless states. To show the 
existence of new orders in gapless states, a new univer- 
sal quantity - projective symmetry group (PSG) - was 
introduced.^" It was argued that (some) PSGs are robust 
against any local perturbations of the Hamiltonian that 
do not change the symmetry of the Hamiltonian. ^^""^^ So 
through PSG, we establish the existence of new orders 
even in gapless states. The new orders are called quan- 
tum order to indicate that the new orders are related 
to patterns of quantum entanglement in the many-body 
ground state. '^^ 



II. SHORT-RANGE AND LONG-RANGE 
QUANTUM ENTANGLEMENT 

What is missed in Landau theory so that it fails to 
describe those new orders? What is the new feature in 
the organization of particles/spins so that the resulting 
order cannot be described by symmetry? 

To answer those questions, let us consider a simple 
quantum system which can be described with Landau 
theory - the transverse field Ising model in two dimen- 
sions: H = —B'Y^Xi — J'Y^ZiZj, where Xi, Yi, and 
Zi are the Pauli matrices on site i. In B ^ J limit, 
the ground state of the system is an equal-weight su- 
perposition of all possible spin-up and spin-down states: 
|$+) = J2{ai} !{'''*})' where {ai} label a particular spin- 
up {Zi = 1) and spin-down {Zi = —1) configuration. In 
the J ^ B limit, the system has two degenerate ground 



states 1$^) = in ■ ■ • t) and |$^) = | || . . . 

The transverse field Ising model has a Z ^ —Z sym- 
metry. The ground state |$^) respect such a symmetry 
while the ground state |$^) (or 1$"^)) break the symme- 
try. Thus the small J ground state |<I>+) and the small 
B ground state 1$^) describe different phases since they 
have different symmetries. 

We note that \^'^) is the exact ground state of the 
transverse field Ising model with B ~ 0. The state has 
no quantum entanglement since j^^) is a direct product 
of local states: 1$''') = t)i where | t)i is an up-spin 
state at site i. The state |$+) is the exact ground state 
of the transverse field Ising model with J = 0. It is also 
a state with no quantum entanglement: |$+) = (g)j(| f 
)i + I oc <»i\+}i where \+)i ee | t)^ + | l)i is a state 
with spin in a;-direction at site i. 

We see that the states (or phases) described by Landau 
symmetry breaking theory has no quantum entanglement 
at least in the J — or B — limits. In the J/B <^ 1 
and B / J <^ 1 limits, the two ground states of the two 
limits still represent two phases with different symme- 
tries. However, in this case, the ground states are not 
unentangled. On the other hand since the ground states 
in the J = or B = limits have finite energy gaps and 
short-range correlations, a small J or a small B can only 
modify the states locally. Thus we expect the ground 
states have only short-range entanglement. 

The above example, if generalized to other symmetry 
breaking states, suggests the following conjecture: // a 
gapped quantum ground state is described by Landau sym- 
metry breaking theory, then it has short-range quantum 
entanglement 7^ 

The direct-product states and short-range entangled 
states only represent a small subset of all possible quan- 
tum many-body states. Thus, according to the point of 
view of the above conjecture, we see that what is missed 
by Landau symmetry breaking theory is long-range quan- 
tum entanglement. It is this long-range quantum entan- 
glement that makes a state to have nontrivial topologi- 
cal/quantum order. 

However, mathematically speaking, the above conjec- 
ture is a null statement since the meaning of short-range 
quantum entanglement is not defined. In the following, 
we will try to find a more precise description (or defi- 
nition) of short-range and long-range quantum entangle- 
ment. We will start with a careful discussion of quantum 
phases and quantum phase transitions. 



III. QUANTUM PHASES AND LOCAL 
UNITARY EVOLUTIONS 

To give a precise definition of quantum phases, let us 
consider a local quantum system whose Hamiltonian has 
a smooth dependence on a parameter g: H{g). The 
ground state average of a local operator O of the sys- 
tem, {0){g), naturally also depend on g. If the function 
{0){g), in the limit of infinite system size, has a singu- 
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larity at gc for some local operators O, then the system 
described by H(g) has a quantum phase transition at 
(7c- After defining phase transition, we can define when 
two quantum ground states belong to the same phase: 
Let |$(0)) be the ground state of H{0) and |$(f )) be the 
ground state of H{1). If we can find a smooth path H{g), 
< g < I that connect the two Hamiltonian H{0) and 
H{1) such that there is no phase transition along the 
path, then the two quantum ground states |<i>(0)) and 
|$(1)) belong to the same phase. We note that "con- 
nected by a smooth path" define an equivalence relation 
between quantum states. A quantum phase is an equiv- 
alence class of such equivalence relation. Such an equiv- 
alence class is called an universality class. 

If |$(0)) is the ground state of H{0) and in the limit 
of infinite system size all excitations above |$(0)) have a 
gap, then for small enough g, we believe that the systems 
described by H{g) are also gapped. '^'^ In this case, we can 
show that, the ground state of H{g), |$((/)), is in the 
same phase as |$(0)), for small enough g, ie the average 
(O) (g) is a smooth function of g near g = for any local 
operator O.'^^ After scaling g to 1, we find that: 

// the energy gap for H{g) is finite for all g in the 
segment [0,1], then there is no phase transition 
along the path g. 

In other words, for gapped system, a quantum phase 
transition can happen only when energy gap closes^^ (see 
Fig. 1)7^ Here, we would like to assume that the reverse 
is also true: a closing of the energy gap for a gapped state 
always induces a phase transition. Or more precisely 



// two gapped states |*&(0)) and |<I'(1)) are in the 
same phase, then we can always find a family of 
Hamiltonian H{g), such that the energy gap for 
H{g) are finite for all g in the segment [0, 1], and 
[^(O)) and |<I'(1)) are ground states of H(Q) and 
H[l) respectively. 



The above two boxed statements imply that two gapped 
quantum states are in the same phase |$(0)) ~ if 
and only if they can be connected by an adiabatic evolu- 
tion that does not close the energy gap. 

Given two states, |<I'(0)) and |$(1)), determining the 
existence of such a gapped adiabatic connection can be 
hard. We would like to have a more operationally practi- 
cal equivalence relation between states in the same phase. 
Here we would like to show that two gapped states 
and |$(1)) are in the same phase, if and only if they are 
related by a local unitary (LU) evolution. We define 
a local unitary(LU) evolution as an unitary operation 
generated by time evolution of a local Hamiltonian for a 
finite time. That is, 

|$(1)) ^ |$(0)) iff |$(1)) =r[e-*^''^9^(9)]|$(0)) (1) 

where T is the path-ordering operator and H{g) = 
TliiOi{g) is a sum of local Hermitian operators. Note 
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FIG. 1: (Color online) Energy spectrum of a gapped sys- 
tem as a function of a parameter s in the Hamiltonian. (a,b) 
For gapped system, a quantum phase transition can happen 
only when energy gap closes, (a) describes a first order quan- 
tum phase transition (caused by level crossing), (b) describes 
a continuous quantum phase transition which has a contin- 
uum of gapless excitations at the transition point, (c) and 
(d) cannot happen for generic states. A gapped system may 
have ground state degeneracy, where the energy splitting be- 
tween the ground states vanishes when system size L — >■ oo: 
VmiL-ioo e = 0. The energy gap A between ground and excited 
states on the other hand remains finite as L — >■ oo. 

that H{g) is in general different from the adiabatic path 
Ii{g) that connects the two states. 

First, assume that two states |$(0)) and |'&(1)) are in 
the same phase, therefore we can find a gapped adia- 
batic path H{g) between the states. The existence of a 
gap prevents the system to be excited to higher energy 
levels and leads to a local unitary evolution, the Quasi- 
adiabatic Continuation as defined in Ref. 34, that maps 
from one state to the other. That is, 

|$(1)) = L/|$(0)), U ^T[e-'iod9H(9)] (2) 

The exact form of H{g) is given in Ref. 33,34 and will be 
discussed in more detail in Appendix. 

On the other hand, the reverse is also true: if two 
gapped states |^'(0)) and |<I'(1)) are related by a local uni- 
tary evolution, then they are in the same phase. Since 
|$(0)) and !$(!)) are related by a local unitary evolution, 
wehave|$(l)) = rie^^-Z'o '*9ff(9)]|$(o)). Let us introduce 

|$(5)) = [/(s)|$(0)), U{s) - T[e-'^S'^9H{a\ (3) 

Assume |<I'(0)) is a ground state of H{0), then |fE'(s)) is 
a ground state of H{s) — U{s)HW (s). If H(s) remains 
local and gapped for all s £ [0, 1], then we have found an 
adiabatic connection between |<i>(0)) and |<I>(1)). 

First, let us show that H{s) is a local Hamiltonian. 
Since if is a local Hamiltonian, it has a form H = J^i Oi 
where Oi only acts on a cluster whose size is ^. ^ is called 
the range of interaction of H. We see that H{s) has a 
form H{s) = J2i Oi{s), where Oi{s) = U{s)OiU\s). To 



show that Oi(s) only acts on a cluster of a finite size, 
we note that for a local system described by H{g), the 
propagation velocities of its excitations have a maximum 
value Vfnax- Since Oi{s) can be viewed as the time evo- 
lution of Oi by H{t) from t — to t — s, we find that 
Oi{s) only acts on a cluster of size ^ + ^ + symax,^"^'^^ 
where ^ is the range of interaction of H. Thus H{s) are 
indeed local Hamiltonian. 

If H has a finite energy gap, then H{s) also have a 
finite energy gap for any s. As s goes for to 1, the 
ground state of the local Hamiltonians, H{s), goes from 
|$(0)) to |$(1)). Thus the two states |$(0)) and |4>(1)) 
belong to the same phase. This completes our argument 
that states related by a local unitary evolution belong to 
the same phase. 

Thus through the above discussion, we show that 



Two gapped ground states,''^ 1^(0)) and |$(1)), 
belong to the same phase if and only if they arc 
related by a local unitary evolution (1). 



A more detailed and more rigorous discussion of this 
equivalence relation is given in Appendix A, where ex- 
act bounds on locality and transformation error is given. 

The relation (1) defines an equivalence relation be- 
tween |$(0)) and The equivalence classes of 
such an equivalence relation represent different quantum 
phases. So the above result implies that the equivalence 
classes of the LU evolutions are the universality classes 
of quantum phases for gapped states. 



IV. TOPOLOGICAL ORDER IS A PATTERN OF 
LONG-RANGE ENTANGLEMENT 
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FIG. 2: (Color online) (a) A graphic representation of a 
quantum circuit, which is formed by (b) unitary operations 
on patches of finite size I. The green shading represents a 
causal structure. 



Since a direct-product state is a state with trivial topo- 
logical order, we see that a state with a short-range en- 
tanglement also has a trivial topological order. This 
leads us to conclude that a non-trivial topological or- 
der is related to long-range entanglement. Since two 
gapped states related by a LU evolution belong to the 
same phase, thus two gapped states related by a local 
unitary evolution have the same topological order. In 
other words. 



Topological order describes the equivalent classes 
defined by local unitary evolutions. 



Or more pictorially, topological order is a pattern of 
long-range entanglement. In Ref. 36, it was shown that 
the "topologically non-trivial" ground states, such as the 
toric code,^^ cannot be changed into a "topologically triv- 
ial" state such as a product state by any unitary locality- 
preserving operator. In other words, those "topologically 
non-trivial" ground states have long-range entanglement. 



Using the LU evolution, we can obtain a more precise 
description (or definition) of short-range entanglement: 



THE LU EVOLUTIONS AND QUANTUM 
CIRCUITS 



A state has only short-range entanglement if and 
only if it can be transformed into an unentangled 
state (ie a direct-product state ) through a local uni- 
tary evolution. 



If a state cannot be transformed into an unentangled 
state through a LU evolution, then the state has long- 
range entanglement. We also see that 



All states with short-range entanglement can 
transform into each other through local unitary 
evolutions. 



Thus all states with short-range entanglement belong to 
the same phase. The local unitary evolutions we con- 
sider here do not have any symmetry. If we require cer- 
tain symmetry of the local unitary evolutions, states with 
short-range entanglement may belong to different sym- 
metry breaking phases, which will be discussed in section 
VI. 



The LU evolutions introduced here is closely related 
to quantum circuits with finite depth. To define quan- 
tum circuits, let us introduce piece-wise local unitary 
operators. A piece-wise local unitary operator has a 
form Upwi = Hi where {Ui\ is a set of unitary op- 
erators that act on non overlapping regions. The size 
of each region is less than some finite number I. The 
unitary operator Upwi defined in this way is called a 
piece- wise local unitary operator with range I. A quan- 
tum circuit with depth M is given by the product of M 
piece- wise local unitary operators (see Fig. 2): t/^^^, = 

^pwi^pwi ' " ^pwi^ ■ quantum information theory, it 
is known that finite time unitary evolution with local 
Hamiltonian (LU evolution defined before) can be sim- 
ulated with constant depth quantum circuit and vice- 
verse. Therefore, the equivalence relation eqn. (f ) can be 
equivalently stated in terms of constant depth quantum 
circuits: 



i$(i)) ^ i$(o)) iff i$(i)) = u^Lim) 



(4) 
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where M is a constant independent of system size. Be- 
cause of their equivalence, we will use the term "Local 
Unitary Transformation" to refer to both local unitary 
evolution and constant depth quantum circuit in general. 

The LU transformation defined through LU evolution 
(1) is more general. It can be easily generalized to study 
topological orders and quantum phases with symmetries 
(see section VI). The quantum circuit has a more 
clear and simple causal structure. However, the quantum 
circuit approach breaks the translation symmetry. So it 
is more suitable for studying quantum phases that do not 
have translation symmetry. 

In fact, people have been using quantum circuits to 
classify many-body quantum states which correspond to 
quantum phases of matter. In Ref. 38, the local uni- 
tary transformations described by quantum circuits was 
used to define a renormalization group transformations 
for states and establish an equivalence relation in which 
states are equivalent if they are connected by a local uni- 
tary transformation. Such an approach was used to clas- 
sify ID matrix product states. In Ref. 39, the local uni- 
tary transformations with disentanglers was used to per- 
form a renormalization group transformations for states, 
which give rise to the multi-scale entanglement renormal- 
ization ansatz (MERA) in one and higher dimensions. 
The disentanglers and the isometrics in MERA can be 
used to study quantum phases and quantum phase tran- 
sitions in one and higher dimensions. Later in this paper, 
we will use the quantum circuit description of LU trans- 
formations to classify 2D topological orders through clas- 
sifying the fixed-point LU transformations. 



VI. SYMMETRY BREAKING ORDERS AND 
SYMMETRY PROTECTED TOPOLOGICAL 
ORDERS 

In the above discussions, we have defined phases with- 
out any symmetry consideration. The H{g) or Upwi in 
the LU transformation does not need to have any sym- 
metry and can be sum/product of any local operators. In 
this case, two Hamiltonians with an adiabatic connection 
are in the same phase even if they may have different sym- 
metries. Also, all states with short-range entanglement 
belong to the same phase (under the LU transformations 
that do not have any symmetry). 

On the other hand, we can consider only Hamiltoni- 
ans H with certain symmetries and define phases as the 
equivalent classes of symmetric local unitary transforma- 
tions: 



I*) 



- i dg W) or I ^ L/f,, I *) 



where H{g) or J7*^c has the same symmetries as H . We 
note that the symmetric local unitary transformation in 

the form T always connect to the identity 

transformation continuously. This may not be the case 
for the transformation in the form U^l^^. To rule out that 
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FIG. 3: (Color online) (a) The possible phases for a Hamil- 
tonian H{gi,g2) without any symmetry, (b) The possible 
phases for a Hamiltonian Hsymmigi, 92) with some symme- 
tries. The shaded regions in (a) and (b) represent the phases 
with short range entanglement {ie those ground states can 
be transformed into a direct product state via a generic LU 
transformations that do not have any symmetry.) 



possibility, we define symmetric local unitary transforma- 
tions as those that connect to the identity transformation 
continuously. 

The equivalent classes of the symmetric LU transfor- 
mations have very different structures compared to those 
of LU transformations without symmetry. Each equiva- 
lent class of the symmetric LU transformations is smaller 
and there are more kinds of classes, in general. 

Fig. 3 compares the structure of phases for a system 
without any symmetry and a system with some symme- 
try in more detail. For a system without any symme- 
try, all the short-range-entangled (SRE) states {ie those 
ground states can be transformed into a direct product 
state via a generic LU transformations that do not have 
any symmetry) are in the same phase (SRE in Fig. 3(a)). 
On the other hand, long range entanglement (LRE) can 
have many different patterns that give rise to different 
topological phases (LRE I and LRE 2 in Fig. 3(a)). 
The different topological orders usually give rise to quasi- 
particles with different fractional statistics and fractional 
charges. 

For a system with some symmetries, the phase struc- 
ture can be much more complicated. The short-range- 
entangled states no longer belong to the same phase, 
since the equivalence relation is described by more spe- 
cial symmetric LU transformations: 

(A) States with short range entanglement belong to dif- 
ferent equivalent classes of the symmetric LU transfor- 
mations if they have different broken symmetries. They 
correspond to the symmetry-breaking (SB) short-range- 
entanglement phases SB-SRE 1 and SB-SRE 2 in Fig. 
3(b). They are Landau's symmetry breaking states. 

(B) States with short range entanglement can belong to 
different equivalent classes of the symmetric LU trans- 
formations even if they do not break any symmetry of 
the system. (In this case, they have the same symme- 
try.) They correspond to the symmetric (SY) short- 
range-entangled phases SY-SRE 1 and SY-SRE 2 in Fig. 
3(b). We say those states have symmetry protected topo- 
logical orders. Haldane phase"^" and Sz = phase of spin- 
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FIG. 4: (Color online) A finite depth quantum circuit can 
transform a state j$) into a direct-product state, if and only if 
the state 1$) has no long-range quantum entanglement. Here, 
a dot represents a site with physical degrees of freedom. A 
vertical line carries an index that label the different physical 
states on a site. The presence of horizontal lines between dots 
represents quantum entanglement. 



1 1 1 



1 1 1 



• 



• 



• 



FIG. 5: (Color online) A piece-wise local unitary transforma- 
tion can transform some degrees of freedom in a state $) into 
a direct product. Removing/adding the degrees of freedom 
in the form of direct product defines an additional equiva- 
lence relation that defines the topological order (or classes of 
long-range entanglement). 



1 chain are examples of states with the same symmetry 
which belong to two different equivalent classes of sym- 
metric LU transformations (with parity symmetry). "^^''^^ 
Band and topological insulators''^"^^ are other examples 
of states that have the same symmetry and at the same 
time belong to two different equivalent classes of symmet- 
ric LU transformations (with time reversal symmetry). 

Also, for a system with some symmetries, the long- 
range-entangled states are divided into more classes 
(more phases): 

(C) Symmetry breaking and long range entanglement can 
appear together in a state, such as SB-LRE 1, SB-LRE 2, 
etc in Fig. 3(b). The topological superconducting states 
are examples of such phases. "'^■^^ 

(D) Long-range-entangled states that do not break any 
symmetry can also belong to different phases such as 
the symmetric long-range-entanglement phases SY-LRE 
1, SY-LRE 2, etc in Fig. 3(b). The many different 
Z2 symmetric spin liquids with spin rotation, transla- 
tion, and time-reversal symmetries are examples of those 
phases. Some time-reversal symmetric topologi- 
cal orders were also called topological Mott-insulators or 
fractionalized topological insulators. ^^"^^ 

VII. LOCAL UNITARY TRANSFORMATION 
AND WAVE FUNCTION RENORMALIZATION 

After defining topological order as the equivalent 
classes of many-body wave functions under LU transfor- 
mations, we like to ask: how to describe (or label) the 
different equivalent classes {ie the different topological 
orders or patterns of long-range entanglement)? 

One simple way to do so is to use the full wave func- 
tion which completely describe the different topological 
orders. But the full wave functions contain a lot of non- 
universal short-range entanglement. As a result, such a 
labeling scheme is a very inefhcient many-to-one labeling 
scheme of topological orders. To find a more efficient or 
even one-to-one labeling scheme, we need to remove the 
non- universal short-range entanglement. 

As the first application of the notion of LU transfor- 
mation, we would like to describe a wave function renor- 
malization group flow introduced in Ref. 58,39. The wave 
function renormalization can remove the short-range en- 



tanglement and simplify the wave function. In Ref. 58, 
the wave function renormalization for string-net states is 
generated by the following two basic moves 

1>( TOT ) ='^*^-* (■■"^^) (5) 

\ ' n ^ ^ 

(Note that the definition of the F-tensor in Ref. 58 is 
slightly different from the definition in this paper.) The 
two basic moves can generate a generic wave function 
renormalization which can reduce the string-net wave 
functions to very simple forms. ^^'^^ Later in Ref. 39, the 
wave function renormalization for generic states was dis- 
cussed in a more general setting, and was called MERA. 
The two basic string-net moves (5) and (6) correspond to 
the isometry and the disentangler in MERA respectively. 
In the MERA approach, the isometrics and the disentan- 
glers are applied in a layered fashion, while in the string- 
net approach, the two basic moves can be applied arbi- 
trarily. In this section, we will follow the MERA setup 
to describe the wave function renormalization. Later in 
this paper, we will follow the string-net setup to study 
the fixed-point wave functions. 

Note that we can use a LU transformation U to trans- 
form some degrees of freedom in a state into direct prod- 
uct (see Fig. 5). We then remove those degrees of free- 
dom in the form of direct product. Such a procedure does 
not change the topological order. The reverse process of 
adding degrees of freedom in the form of direct prod- 
uct also does not change the topological order. We call 
the local transformation in Fig. 5 that changes the de- 
grees of freedom a generalized local unitary (gLU) trans- 
formation. It is clear that a generalized local unitary 
transformation inside a region A does not change the re- 
duced density matrix pA for the region A. This is the 
reason why we say that (generalized) local unitary trans- 
formations cannot change long-range entanglement and 
topological order. Similarly, the addition or removal of 
decoupled degrees of freedom to or from the Hamilto- 
nian, H ^ H (E) H^^^ will not change the phase of the 
Hamiltonian (ie the ground states of H and H ® Hdp are 



8 



1 1 i 1 i 



region A 



(a) 



VIII. WAVE FUNCTION RENORMALIZATION 
AND A CLASSIFICATION OF TOPOLOGICAL 
ORDER 



\ 7 
I I I 

i 11 i 



region A 



(b) 



FIG. 6: (Color online) (a) A gLU transformation U acts in 
region A of a state |$), which reduces the degree freedom in 
region A to those contained only in the support space of 1$) 
in region A. (b) U'^U — P is a. projector that does not change 
the state 1$). 



in the same phase) , if those degrees of freedom form a di- 
rect product state (ie the ground state of Hdp is a direct 
product state). 

Let us define the gLU transformation U more carefully 
and in a more general setting. Consider a state |$). Let 
Pa be the reduced density matrix of |$) in region A. pA 
may act in a subspace of the total Hilbert space Va in 
region A, which is called the support space of region 
A. The dimension D'^J of V^' is called support dimension 
of region A. Now the Hilbert space Va in region A can be 
written as Va — V^ ® ^a ■ Let \ipi), « = 1, ■■■,0"/ be a 
basis of this support space V^, \tpi), i = D'^/ + 1, Da 
be a basis of V^^, where Da is the dimension of Va, and 
i = 1,...,Da be a basis of Va- We can introduce 
a LU transformation [Z-^"" which rotates the basis Itpi) 
to \ipi). We note that in the new basis, the wave func- 
tion only has non-zero amplitudes on the first basis 
vectors. Thus, in the new basis \ipi), we can reduce the 
range of the label i from [1, Da] to [1, D^'] without losing 
any information. This motivates us to introduce the gLU 
transformation as a rotation from \ipi), i = ^,---,Da to 
\ipi), « = 1, D'^A ■ "L^^ rectangular matrix U is given by 
Uij — {ipi\^j). We also regard the inverse of W, as 
a gLU transformation. A LU transformation is viewed 
as a special case of gLU transformation where the de- 
grees of freedom are not changed. Clearly U^U = P and 
UW = P' are two projectors. The action of P does not 
change the state |$) (see Fig. 6(b)). 

We note that despite the reduction of the degrees of 
freedom, a gLU transformation defines an equivalent rela- 
tion. Two states related by a gLU transformation belong 
to the same phase. The renormalization flow induced by 
the gLU transformations always flows within the same 
phase. 



As an application of the wave function renormalization, 
in this section, we will study the structure of fixed-point 
wave functions under the wave function renormalization, 
which will lead to a classification of topological order 
(without any symmetry). 

We note that as wave functions flow to a fixed point, 
the gLU transformations in each step of the renormal- 
ization also flow to a fixed point. So instead of study- 
ing fixed-point wave functions, here, we will study the 
fixed-point gLU transformations. For this purpose, we 
need to fix the renormalization scheme. In the following, 
we will discuss a renormalization scheme motivated by 
the string-net wave function. After we specify a proper 
wave function renormalization scheme, then the fixed- 
point wave function is simply the wave function whose 
"form" does not change under the wave function renor- 
malization. 

Since those fixed-point gLU transformations do not 
change the fixed-point wave function, their actions on the 
fixed-point wave function do not depend on the order of 
the actions. This allows us to obtain many conditions 
that gLU transformations must satisfy. From those con- 
ditions, we can determine the forms of allowed fixed-point 
gLU transformations. This leads to a general description 
and a classification of topological orders and their corre- 
sponding fixed-point wave functions. 

The renormalization scheme that we will discuss was 
first used in Ref. 58 to characterize the scale invariant 
string-net wave function. It is also used in Ref. 17 to sim- 
plify the string-net state in a region, which allows us to 
calculate the entanglement entropy of the string-net state 
exactly. A similar approach was used in Ref. 59 to show 
quantum-double/string- net wave function to be a fixed- 
point wave function and its connection to 2D MERA.'^^ In 
the following, we will generalize those discussions by not 
starting with string-net wave functions. We just try to 
construct local unitary transformations at a fixed point. 
We will see that the fixed-point conditions on the gLU 
transformations lead to a mathematical structure that is 
similar to the tensor category theory - the mathematical 
framework behind the string-net states. 



A. Quantum states on a graph 

Since the wave function renormalization may change 
the lattice structure, we will consider quantum states de- 
fined on a generic trivalence graph G: Each edge has 
iV + 1 states, labeled by i = 0,...,N (see Fig. 7). We 
assume that the index i on the edge admits an one-to-one 
mapping *: i ^ i* that satisfies (i*)* = i. As a result, 
the edges of the graph are oriented. The mapping i i* 
corresponds to the reverse of the orientation of the edge 
(see Fig. 8). Each vertex also has physical states, labeled 
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FIG. 7: A quantum state on a graph G. There are A''+l states 
on each edge which are labeled by / = 0, A*'. There are A'^„ 
states on each vertex which are labeled by a = 1, A'^i,. 




FIG. 9: A quantum state on a graph G. If the three edges 
of a vertex are in the states i, j, and k respectively, then the 
vertex has Nijk states, labeled by a = 1, Nijk- Note the 
orientation of the edges are point towards to vertex. Also 
note that i j ^ k runs anti-clockwise. 



FIG. 8: The mapping i ^ i* corresponds to the reverse of 
the orientation of the edge 



by a = l,...,iV„ (see Fig. 7). 

Each labeled graph (see Fig. 7) corresponds to a state 
and all the labeled graphs form an orthonormal basis. 
Our fixed-point state is a superposition of those basis 



states: |$fix) = J2 



all . 



fix 



) 




Here we will make an important assumption about the 
fixed-point wave function. We will assume that the fixed- 
point wave function is "topological" : two labeled graphs 
have the same amplitude if the two labeled graphs can 
be deformed into each other continuously on the plane 
without the vertices crossing the links. For example 

TT^^X / ^ 





V'fix I " '^^"^ 

such an assumption, the topological orders studied in this 
paper may not be most general. 

We also assume that all the fixed-point states on each 
different graphs to have the same "form" . This assump- 
tion is motivated by the fact that during wave function 
renormalization, we transform a state on one graph to a 
state on a different graph. The "fixed-point" means that 
the wave functions on those different graphs are all deter- 
mined by the same collection of the rules, which defines 
the meaning of having the same "form" . However, the 
wave function for a given graph can have different total 
phases if the wave function is calculated by applying the 
rules in different orders. Those rules are noting but the 
fixed-point gLU transformations. 



B. The structure of entanglement in a fixed-point 
wave function 



Before describing the wave function renormalization, 
let us examine the structure of entanglement of a fixed- 
point wave function. First, let us consider a fixed-point 
wave function <f>fix on a graph. We examine the wave 



function on a patch of the graph, for example. 



The fixed-point wave function 




(only the 



relevant part of the graph is drawn) can be viewed as a 
function of a,/?, to: Viifc;,r(a, ™) = *fix ^^i^j^^ if 

we fix I, J, fc, I and the indices on other part of the graph. 
(Here the indices on other part of the graph is summa- 
rized by F.) As we vary the indices F on other part of 
graph (still keeping i, k, and I fixed), the wave func- 
tion of a,/3,TO, ij^iju^vioi, j3,m), may change. All those 
'0ijfci,r(«) /?) w) form a linear space of dimension Diji^i*. 
Dijki is an important concept that will appear later. We 
note that the two vertices a and /? and the link m form 
a region surrounded by the links k, I. So we will call 
Due to the dimension-Di^/;.;. space the support space Vijki* and 



Dijki' the support dimension for the state $fix on the 
region surrounded by the fixed boundary state i,j,k,l. 
Similarly, we can define Dijk as the support dimension 



of the $fix 




on a region bounded by links i, j, k. 



Since the region contains only a single vertex a with 
physics states, we have Dij}^ < iV„. We can use a lo- 
cal unitary transformation on the vertex to reduce the 
range of a to 1, ...,Nijk where Nijk = Dijk- In the rest 
of this paper, we will implement such a reduction. So, 
the number of physical states on a vertex depend on the 
physical states of the edges that connect to the vertex. 
If the three edges of a vertex are in the states i, j, and k 
respectively, then the vertex has Nijk states, labeled by 



a 



1, 



ijk 



(see Fig. 9). Here we assume that 



N,jk = N 



jki- 



(7) 



We note that in the fixed-point wave function 
'I'fix I I' the number of choices of a,/3,m is 




Nijki' = J2m=o ^jim' ^kmi' ■ Thus the support dimen- 
sion Dijki* satisfies Dijki* < Nijki*. Here we will make 
an important assumption - the saturation assumption: 
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The fixed-point wave function saturate the inequality: 

N 



(8) 



We will see that the entanglement structure described by 
such a saturation assumption is invariant under the wave 
function renormalization. 



where ~ means equal up to a constant phase factor. 
(Note that the total phase of the wave function is unphys- 
ical.) We will call such a wave function renormalization 
step a F-move. 

We would like to remark that eqn. (11) relates two wave 
functions on two graphs Gi and G2 which only differ by a 
local reconnection. We can choose the phase of F^-'^^'^f 
to make ~ into —: 



C. The first type of wave function renormalization 

Our wave function renormalization scheme contains 
two types of renormalization. The first type of renormal- 
ization does not change the degrees of freedom and corre- 
sponds to a local unitary transformation. It corresponds 

to locally deform the graph ^^^^ to . (The parts 

that are not drawn are the same.) The fixed-point wave 
function on the new graph is given by $fix \§'^ 
Again, such a wave function can be viewed as a function 
ofx,S,n: '4'ijki,rix,S,n) = $fix I \^ I if we fix j, fc, ? 
and the indices on other part of the graph. The support 
dimension of the state <i>fix on the region sur- 

rounded by i,j,kj is D^jki'-- Again Ajfc/* < N.jki*, 
where Nijki* = Y,n=o ^kjn'Ni.m is number of choices of 
X, 6, n. The saturation assumption implies that Nijki* = 
Dijki* ■ 

The two fixed-point wave functions <i>fix ^^*m^^^^ ^^'^ 

'i'fix are related via a local unitary transforma- 

tion. Thus Dijki* = Dijki" , which implies 

N N 
J2 NJ^„^* Nkral' = Nkjn' ^i'm • (9) 



m— n— 

We express such an unitary transformation as 

n=0 x=l S=l 

(10) 

or graphically as 

''i j k\ N N„,i. (i j y k\ 

i>fix I ^ 1 - E E E ^^LT/'f- 




n=0 x=l <5=1 



(11) 





E E E 

n=0 x=l <5=1 



But such choice of phase only works for a particular pair 
of graphs Gi and G2. To use Flj^l^f to relate all pair of 
states that only differ by a local reconnection, in general, 
we may have a phase ambiguity, with the value of phase 
depend on the pair of graphs. So eqn. (11) can only be 
a relation up to a total phase factor. 

For fixed i, j, k, and I, the matrix i*"^; with ma- 
trix elements {FlD'^^^g^ = Fll'^^'^f is a N.^ki- dimen- 
sional matrix (see (9)). The mapping '0ijfe;.r(Xi '^i ~^ 



'^ijki.ri.Oi, I3,m) generated by the matrix F^j is unitary. 
Since, as we change F, rpijki,rix, <5, n) and ipijki.r{a, (3, m) 
span two Nijki' dimensional spaces. Thus F^j is a 



N, 



■ijkl' 



X ^ijki'^ unitary matrix 



E 



F. 



ma/B ,m' a' 13' i 



(12) 



where 5mai3,m'a'p' ~ 1 when m = to', a — a' , /S — /?', 
and S-fnap,m'a'i3' = Otherwise. 



We can deform 
We see that 



to '^<<^, and \^ to 



N 



^-1 V i-EE^'X.r'j'-iN^^ 

V / n=0 x,S 



(13) 



Eqn. (11) and eqn. (13) relate the same pair of graphs, 
and thus 



rpijm.ap ^ Tjikl* m* .13a. 
kln,xS ij*n* ,Sx 



(14) 



(where we have used the condition Nijki = Dijki-) 
Using the relation (12), we can rewrite eqn. (11) as 



<i>fix I \j>^ - E E(^^LT/)*'J'fi>^ 

'/ / m=0 a,l3 



(15) 



11 



We can also express $fix as 




(16) 



When Njim' < 1 or Nkmi' < 1, the left-hand-side of 
eqn. (11) is always zero. Thus F^j^^'^f = when Njim* < 
1 or Nkmi' < 1- When Nkjn' < 1 or N^ii* < 1, wave 
function on the right-hand-side of eqn. (11) is always 



zero. So we can choose F, 



or Nr,ii* < 1. 



when N^jn' < 1 



using the relabeled eqn. (11). So we have 

\^kln,xS > — ^l'i'm',l3a- \^ ' I 

Since the total phase of the wave function is unphysical, 
the total phase of F^jl^'^f can be chosen arbitrarily. We 

can choose the total phase of F^j^:^/ to make 

y^kln,xS > ~ ^I'i'm'^Pa- > 

If we apply eqn. (18) twice, we reproduce eqn. (14). Thus 
eqn. (14) is not independent and can be dropped. 
From the graphic representation (11), We note that 



^:LT/ = -hen 



(19) 



Njira' < 1 or Nkrnl' < 1 Or Nkjn' < 1 Or Nml' < 1. 



The F-move (11) maps the wave functions on two dif- 
ferent graphs through a local unitary transformation. 
Since we can locally transform one graph to another 
graph through different paths, the F-move (11) must sat- 
isfy certain self consistent condition. For example the 



graph 




can be transformed to 




through 



p p 
two different paths; one contains two steps of local trans- 
formations and the other contains three steps of local 
transformations as described by eqn. (11). The two paths 
lead to the following relations between the wave func- 
tions: 




ETpinkn,l3xPh 



V F 



qps,4>-f 




(20) 



fix 




V ^ / 



rjip 



V p J 





T^ijm.Oi^ T^itn,<^X Tpj}^t,r]K.^ 



^ ^ ^ knt^rjif ^ Ips.K'j ^ Isq 




V " J 



(21) 



The consistence of the above two relations leads a condi- 
tion on the F tensor. 

To obtain such a condition, let us fix i, j, fc, I, p. 




and view <i>fix 



■4){a,(3,x,m,n) = $fix 



as a function of a, /3, x, "ni, n: 
As we vary indices 




on other part of graph, we obtain different wave func- 
tions "0(0;, /?, X, m, n) which form a dimension Dijkip- 
space. In other words, Dijkip* is the support dimension 
of the state $fix on the region a, /?, x, rn, n with bound- 
ary state i,j,k,l,p fixed (see the discussion in section 
VllIB). Since the number of choices of a, /3,x,™,"' is 

N^jklp* = J2m,n ^jim' Nkmn' Ninp' , WC havC D^jklp' < 

Nijkip' ■ Here we require a similar saturation condition 
as in (8): 



ijklp* 



D 



ijklp* 



(22) 
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Similarly, the number of choices of S, (f), 7, q, s in 



/'■ i k i\ 




is also N, 



ijklp* ■ 



Here we again assume 



• , where Dijkip* is the support dimension 



on the region bounded by i, j, fc, I, p. 



So the two relations (20) and (21) can be viewed as 
two relations between a pair of vectors in the two -Dijfcip* 
dimensional vector spaces. As we vary indices on other 
part of graph (still keeping i,j,k,l,p fixed), each vector 
in the pair can span the full Dij^ip- dimensional vector 
space. So the validity of the two relations (20) and (21) 
implies that 

^fcjt* Ntir.* Nu,. 

/ J / J / J knt,7]ip Ips.K^ lsq.54> 
t rj—1 K— 1 



A6p ST^ p'rnkn,l3x pijm,a 
/ J lpq,Se qps,<t>J 



(23) 




.0 




k ^ \k 



FIG. 10: A "triangle" graph can be transformed into a "tad- 
pole" via two steps of the first type of wave function renor- 
malization (ie two steps of local unitary transformations). 



via the first type of renormalization steps (see Fig. 10), 
which are local unitary transformations. In the simplified 



second type renormalization, we want to reduce 



to , so that we still have a trivalence graph. This re- 
quires that the support dimension Dw of the fixed-point 
/ , 

is given by 



wave function $ 




(25) 



which is a generalization of the famous pentagon iden- 
tity (due to the extra constant phase factor e'^'^). We 
will call such a relation projective pentagon identity. The 
projective pentagon identity is a set of nonlinear equa- 
tions satisfied by the rank- 10 tensor F^j^'^f and 9p. The 
above consistency relation is equivalent to the require- 
ment that the local unitary transformations described by 
eqn. (11) on different paths all commute with each other 
up to a total phase factor. 



D. The second type of wave function 
renormalization 

The second type of wave function renormalization does 
change the degrees of freedom and corresponds to a gen- 
eralized local unitary transformation. One way to im- 
plement the second type renormalization is to reduce 



(the part of the graph that is not 




drawn is unchanged): 
/ 





(24) 



But we can define a simpler second type renormaliza- 

p >■ 



tion, by noting that 




can be reduced to 
•J i y*- 



This implies that 




(26) 



The simplified second type renormalization can now be 
written as (since Da- = 1) 




(27) 



We will call such a wave function renormalization step a 
P-move.^^ Here p'^^'"^ satisfies 



k,j a=l l3=l 



(28) 



and 



pkj,aP ^ 0^ j^^^^^ ^ ^ ^^^^^^ ^ (29) 



The condition (28) ensures that the two wave functions 
on the two sides of eqn. (27) have the same normalization. 
We note that the number of choices for the four indices 
(j, k, a, P) in p^^'"^ must be equal or greater than 1: 



a 



(30) 
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Notice that 



fix 



m,A,7 




Using eqn. (27) and its variation 



i ii i 



(31) 



fix 



Pi"'^''$fix 



(-) 



(32) 



we can rewrite eqn. (31) as 

^\Opipkj,ap _ p3fk,pOL pi*mj,X^ plm.jiu 

i / J i*i*m* ,X'-y m*-i*l.i'fi i* 

m,A,7.Z,i/./i 

(33) 

which is a condition on p^^'"^ , 

More conditions on F^-'^^l^f and P^'^'°'^ can be ob- 



N 




n=0 x,<5 




(34) 



which imphes that 



i^kiri^x ^ Pi'''"'^)^ in this section, we are going to discuss 
how to calculate the corresponding fixed-point wave func- 
tion $fix from the solved fixed-point gLU transformation 

f-pijm,aP pkj,al3\ 

First we note that, using the two types of wave func- 
tion renormalization introduced above, we can reduce 



any graph to So, once we know <i>fix |^V^/^J > '^6 

can reconstruct the full fixed-point wave function $fix on 
any connected graph. 
Let us assume that 



fix V / 



= A' = A' 



Here A^ satisfy 

A'^A'\ "^A'{A')* = 1. 



(37) 



(38) 



The condition ^^A*(A*)* = 1 is simply the normal- 
ization condition of the wave function. The condition 



A^ = A^ come from the fact that the graph can be 

deformed into the graph (^T) on & sphere. 

To find the conditions that determine A', let us first 
consider the fixed-point wave function where the index on 



a link is i: $fix(i, T) — $fix 



where F are indices 



on other part of graph. We note that the graph 

can be deformed into the graph ^ on a sphere. Thus 

$flx ( ^ — ^ I = 'J' fix ( ""V^ I • Using the F-moves and 



the P-moves, we can reduce 



to 



^klnaS Mc- ''fcn'Pfix 



'i k\ 

\5y 



(35) 



<ffix( 



(39) 



We find 



^ _ T^ijm,al3 pjp,XV 

■5 — ^klk,x& ^k' 

X=l 

for all fc, i, I satisfying Nku* > 0. (36) 



E. The fixed-point wave functions from the 
fixed-point gLU transformations 



In the last section, we discussed the conditions that 
a fixed-point gLU transformation {Fli^^x , P^-^'-'^^) must 
satisfy. After finding a fixed-point gLU transformation 



We see that 



$fix (Q) = 



(40) 



for all i. Otherwise, any wave function with i-link will 
be zero. 

To find more conditions on A^, we note that 




(41) 
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By rotating 



Pj. ' ' A-) or 



by 180°, we can show that p^i'^^A^ 



We also note that 



fix 



. i I m,A,7 




This ahows us to show 



ikj,afi — ^ ^ j'i'm,X^^i ^ 

m, A,7 



(42) 



(43) 



(44) 



where 




(45) 



^fcji ap comes from the fact that 



The condition ^f^^- „^ 

the graph "(^^^^P and the graph "(^^^P can be deformed 



into each other on a sphere. 

Also, for any given k,a that satisfy N]^ji > 0, the 

wave function $fix ^"^^ ^ ^ must be non-zero for some 

r, where T represents indices on other part of the graph. 
Then after some F-moves and P-moves, we can reduce 

"C^'^- foi' s-iiy given i,j,k,a that sat- 
isfy Nkji > 0, $fix ^"(^^^P^ is non-zero for some /?. 

Since such a statement is true for any choices of ba- 
sis on the vertex a, we find that for any given i,j,k^a 
that satisfy N^ji > and for any non-zero vector Va, 

X]q ^ct'J'fix ^"(^^^P^ is non-zero for some /3. This means 
that the matrix Mkji is invertible, where Mkji is a matrix 



The above also implies that 



(47) 



The conditions eqns. (38, 42, 44, 45) allow us to deter- 
mine A' (and „^). 

From eqn. (45), we see that relation $fix ^"(^^^^^P^ — 

$fix ^"(^^^P^ leads to some equations for F^in '^x ■• 

pkj,a^ ^ and A'-. More equations for F^^^'^x ^ pkj,ai3 ^ 
A^, can be obtained by using the relations 




(48) 




from the tetrahedron rotation symmetry ■ and 



7A 



jX.pcre y 



^ T^ijm' ,al3 -pln'i' ,SX ppl* ^ae^ 

— ^kln.jX ^nlp',ae ^fix 

7A,pcre 

Epijm' ,af) pln'i' ,SX ppl* ,a<i^9 
kln,jX nip* ,ae n ^kjn* 

yX^pcre 




(49) 



It is not clear if eqn. (48) and eqn. (49) will lead to new 
independent equations or not. In the following discus- 
sions, we will not include eqn. (48) and eqn. (49). We 
find that, at least for simple cases, the equations without 
eqn. (48) and eqn. (49) are enough to completely deter- 
mine the solutions. 

To summarize, the conditions (9, 30, 47, 12, 18, 23, 33, 
36, 40, 38, 42, 44, 45, 46) form a set of non-linear equa- 



whose elements are given by iMkJ^)a(^ = *fix ■ tions whose variables are N,,k, -Ffef";" , ^ > and 



and A' 



Let us define det 
that 



fix 



aP-HP 



= det(Mfcji), we find 



det 




(46) 



ydF,dpi,8p2)- Finding Nijk, F^^^.^^ , ^ ■ 
that satisfy such a set of non-linear equations corresponds 
to finding a fixed-point gLU transformation that has a 
non-trivial fixed-point wave function. So the solutions 

{Nijk,Fl-i^^x ^Fi'''"^ j^^) give us a characterization of 
topological orders. This may lead to a classification of 
topological order from the local unitary transformation 
point of view. 
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IX. SIMPLE SOLUTIONS OF THE 
FIXED-POINT CONDITIONS 

In this section, let us find some simple solutions of the 
fixed-point conditions (9, 30, 47, 12, 18, 23, 33, 36, 40, 38, 
42, 44, 45, 46) for the fixed-point gLU transformations 
(iVijfe, F^j^;^"'^, P-*^^'"^) and the fixed-point wave function 



A. Unimportant phase factors in the solutions 



tion also does not affect the fusion rule Nij^. The new lo- 
cal unitary transformation changes (Flil^^x ^ ^i'' "^ ^ 
and generates a continuous family of the fixed-point wave 
functions parameterized by the unitary matrix u^^°i°'^°\ 
Again, those fixed-point wave functions all belong to the 
same phase. 

In the following, we will study some simple solutions of 
the fixed-point conditions. We find that, for those exam- 
ples, the solutions have no addition continuous parame- 
ter apart from the phase factors discussed above. This 
suggests that the solutions of the fixed-point conditions 
correspond to isolated zero-temperature phases. 



Formally, the solutions of the fixed-point conditions 
are not isolated. They are parameterized by several con- 
tinuous phase factors. In this section, we will discuss the 
origin of those phase factors. We will see that those dif- 
ferent phase factors do not correspond to different states 
of matter (ie different equivalence classes of gLU trans- 
formations) . So after removing those unimportant phase 
factors, the solutions of the fixed-point conditions are 
isolated (at least for the simple examples studied here). 

We notice that, apart from two normalization condi- 
tions, all of the fixed-point conditions are linear in p^^'"^ 
and A\ Thus ff {Fjfin^:^f , P^'''"'^ , A') is a solution, then 
{KZ:,f , e'^' Pt''''^ ,e''^^A') is also a solution. How- 
ever, the two phase factors e"^^-^ do not lead to different 
fixed-point wave functions, since they only affect the to- 
tal phase of the wave function and are unphysical. Thus 
the total phases of P''-' and A"^ can be adjusted. We can 



,004 1 



> and 



use this degree of freedom to set, say, Pq 
A°>0. 

Similarly the total phase of F^^^'^x unphysical 
and can be adjusted. We have used this degree of free- 
dom to reduce eqn. (17) to eqn. (18). But this does not 

The transforma- 



totally fix the total phase of F, 



tion P, 



kln^'yX 



does not affect eqn. (18). We 
can use such a transformation to set the real part of a 
non-zero component of F^j^^-f to be positive. 

The above three phase factors are unphysical. How- 
ever, the fixed-point solutions may also contain phase 
factors that do correspond to different fixed-point wave 
functions. For example, the local unitary transfor- 
mation e'^'o^^'o does not affect the fusion rule Nij^, 
where Mig is the number of links with |^o)-state and 
|/Q)-state. Such a local unitary transformation changes 

i^kin^^x ' '^i'''"'^ ^ ^^'^ generates a continuous family 
of the fixed-point wave functions parameterized by 9ig. 
Those wave functions are related by local unitary trans- 
formations that continuously connect to identity. Thus, 
those fixed-point wave functions all belong to the same 
phase. 

Similarly, we can consider the following local unitary 
transformation \a) — > J2i3^ap''"''"^\l^) ^^^^ '^'^ each 
vertex with states |io), \ jo), \ko) on the three edges con- 
necting to the vertex. Such a local unitary transforma- 



B. iV = 1 loop state 

Let us first consider a system where there are only two 
states |0) and |1) on each link of the graph. We choose 
i* = i and the simplest fusion rule that satisfies eqn. (9), 
eqn. (30), eqn. (47) is 



Nooo = iVno = ^101 = iVou = 1, 
other Nijk = 0. 



(50) 



Since Nijk < 1, there is no states on the vertices. So 
the indices a, /3, ... labeling the states on a vertex can be 
suppressed. 

The above fusion rule corresponds to the fusion rule 
for the = 1 loop state discussed in Ref. 58. So we will 
call the corresponding graphic state iV = 1 loop state. 

Due to the relation (18), the different components of 
the tensor P^j™ are not independent. There are only four 
independent potentially non-zero components which are 
denoted as fo,...J^: 



p; 



000 , 



pOOO 

^111, 



^011/ 

pllO 
^110 



.1 -io 



100^ 



77.110 

-^001 



> 



= {Fi^i 

= h 



Xr = h 



(51) 



We note that P^^™ 



in eqn. (11) relates wave functions 
on two graphs. In the above we have drawn the two 
related graphs after the P tensor, where the first graph 
following P corresponds to the graph on the left-hand 
side of eqn. (11) and the second graph corresponds to 
the graph on the right-hand side of eqn. (11). The doted 
line corresponds to the |0)-state on the link and the solid 
line corresponds to the |l)-state on the link. There are 
four potentially non-zero components in P^\ which are 
denoted by pov,P3: 



pOO 

-^0 



Po, 



pOl 



pi, pr=P2, Pi"^=P3. (52) 
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We can adjust the total phases of pi and to make 
Po > Q and A° > 0. We can also use the local unitary 
transformation e'^'o*^'o with = 1 to make /i > 0, since 
the i^'s described by /i in eqn. (51) are the only F^s that 
change the number of |1) -links. 

The fixed-point conditions (9, 30, 47, 12, 18, 23, 33, 
36, 40, 38, 42, 44, 45, 46) form a set of non-linear equa- 
tions on the ten variables fi, pi, and A^. Many of the 
non- linear equations are dependent or even equivalent. 
Using a computer algebraic system, we simplify the set 
of non-linear equations. The simplified equations are (af- 
ter making the phase choice described above) 



/2 = 1, /s = V, 
VPo, P2=P0, 



fo — fi 
Pi =P3 

Po + bil' = i- b2|' + b3l' = l, 
p,A°=p,A\ 7^p,A'^p,A', 



(53) 



\A'\' = 1 



where rj — ±1. The above simplified equations can be 
solved exactly. We find two solutions parameterized by 
rj = ±1: 



/o = fi = /2 = 1, /s = V, 
1 

Pa=P2 = Pi=P3 = 



A" 



1 

71 



V2' 
A' 



V 

V2' 



V 



We also find 



1. 



(54) 



(55) 



The T] = 1 fixed-point state corresponds to the Z2 loop 
condensed state whose low energy effective field theory is 
the Z2 gauge theory.^*'^° We call such a state, simply, the 
Z2 state. The ij = —1 fixed-point state corresponds to 
the double-semion state whose low energy effective field 
theory is the C/(l) x [/(I) Chern-Simons gauge theory^^'^° 



(56) 



C. A'^ = 1 string-net state 

To obtain another class of simple solutions, we modify 
the fusion rule to 

A^ooo = Nno = Nwi = A^oii = A^iii = 1, 

other N,,k = 0. (57) 

while keeping everything the same. The above Nijk also 
satisfies eqn. (9), eqn. (30), and eqn. (47). 

The new fusion rule corresponds to the fusion rule for 
the A'' = 1 string-net state discussed in Ref. 58. So we 
will call the corresponding graphic state = 1 string-net 
state. 



Again, due to the relation (18), the different compo- 
nents of the tensor F^fl^ are not independent. Now there 
are seven independent potentially non-zero components 
which are denoted as /o,...,/6: 



^000 ■ 



-^111 



^011 / 
^111 



pllO 
-^110 



^111 

Fill 
^111 



^ l^ioo 

— piio 

— -^001 

— I^IOI 

— l^iii 

— ('pill 

= h 

= /6 



Xr = h 

- Fill 
^ ^011 

= h 
y = h 



(58) 



Note that F's described by /i and arc the only F's 
that change the number of |l)-links and the number of 
|l)|l)|l)-vertices. So we can use the local unitary trans- 
formation e'(^^^i+<^*ii_i) to make /i and /s to be positive 
real numbers. (Here Mi is the total number of |l)-links 
and Mill is the total number of |l)|l)|l)-vertices.) We 
also use the freedom of adjusting the total sign of i^^;^ 
to make Re(/o) > 0. 

There are five potentially non-zero components in P^'' , 
which are denoted by po ■ iP4 ^ 



pOl 



Pa, 
Pa, 



Pi' 



Pi, 
Pa- 



P2, 



(59) 



We use the freedom of adjusting the total phase of P^^ 
to make Pq to be a positive number. We can also use the 
freedom of adjusting the total phase of A^ to make A^ to 
be a positive number. 

The fixed-point conditions (12, 18, 23, 33, 36, 40, 38, 
42, 44, 45, 46) form a set of non-linear equation on the 
variables /i, pi, and A', which can be simplified. The 
simplified equations have the following form 



/o — /i — /2 — /3 — 1, h — /s 



-h > 0, 



Pifl+Pi = 1, Po = hPi, P2 = Pa, P3 = Pi, Pi = 



A^ 



UA\ [A^ + iAy^l, ft + U-1^0 



(60) 



Let 7 be the positive solution of 7^^ + 7 = 1: 7 
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We see that /a = The above can be written as 

/o = /i = /2 = /a = 1, .A = -/e = 7, /s = V7, 



7 1 . 

P0=P2= 2 , Pi = P3 = 2 , -, , P4 = 0, 



7^ + 1 



A' = 



7^ + 1 



(61) 



72 + 1' 
We also find 

The fixed-point state corresponds to the = 1 
string-net condensed state^^ whose low energy effective 
field theory is the doubled SO{3) Chern-Simons gauge 
theory.s° 



D. An N — 2 string-net state — the Z3 state 

The above simple examples correspond to non- 
orientable string-net states. Here we will give an example 
of orientable string-net state. We choose = 2, 0* = 0, 
1* = 2, 2* = 1, and 



There are nine potentially non-zero components in P^'' , 
which are denoted hy po,...,ps- 

PS°=Po, P^'^Pu P^^=P2, Pr=P3, P'.'^PA, 
Pr=P^, PT^P^. Pi'-Pi, P^'-Ps- (65) 

Using the transformations discussed in section IX B, we 
can fix the phases of /i, /2, fe, and po to make them 
positive. 

The fixed-point conditions (12, 18, 23, 33, 36, 40, 38, 
42, 44, 45, 46) form a set of non-linear equation on the 
variables fi, pi, and A^, which can be solved exactly. Af- 
ter fixing the phases using the transformations discussed 
in section IX B, we find only one solution 

/, = 1, z = 0,l,...,7, 
1 



, . . . , 



A" = A' = A^ = 



1 



(66) 



We also find 



Nqoq — Nqi2 — A^120 — ^^201 — ^021 — -^102 — ^210 

= ^111 = iV222 = 1. (63) 

The above iVy ^ satisfies eqn. (9), eqn. (30), and eqn. (47). 

Due to the relation (18), the different components of 
the tensor i^^;™ are not independent. There are eight 
independent potentially non-zero components which are 
denoted as /o,...,//: 



1. 



-^000 , 



^111 ■ 



000 



-f^oii,- - 



fo 

(F 



F, 



120 V 



002/ 



202 



J* - /l 



/p022W 

l-'^ioo ■ 



-^001 



h 



{fxX. 



{Fl^l'^ 



{Fw2 
[^212 



( (7^221 



p221 
^120 

( ^210- 
{^220, 



Xr 
Xr 
Xr 
Xr 
Xr 



h 

h 

iF?i^H Xr 



r = f7 



(64) 



(67) 



The fixed-point state corresponds to the Z3 string-net 
condensed state"''^ whose low energy effective field theory 
is the U{1) X U{1) Chern-Simons gauge theory^^''''* 



C 



(68) 



which is the Z3 gauge theory. 

We note that all the above simple solutions also satisfy 
the standard pentagon identity, although we solved the 
weaker projective pentagon identity. It is not clear if 
we can find non-trivial solutions that do not satisfy the 
standard pentagon identity. 



X. A CLASSIFICATION OF TIME REVERSAL 
INVARIANT TOPOLOGICAL ORDERS 

There are several ways to define time reversal operation 
for the graphic states. The simplest one is given by 



T : $(r) $*(r) 



(69) 



where F represents the labels on the vertices and links 
which are not changed under T. (This corresponds to 
the situation where the different states on the links and 
the vertices are realized by different occupations of scalar 
bosons.) For such a time reversal transformation, — 1 
and the real solutions of the fixed-point conditions (9, 
30, 47, 12, 18, 23, 33, 36, 40, 38, 42, 44, 45, 46) give 
us a classification of time reversal invariant topological 
orders in local spin systems. Note that the time reversal 
invariant topological orders are equivalent class of local 
orthogonal transformations that connect to the identity 
transformation continuously. 
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Different real solutions (A^i.fe, -F^QT, i^"''"'', of 
the fixed-point conditions do not always correspond 
to different time reversal invariant topological orders. 
The solutions differ by some unimportant phase factors 
(which are ±1 signs) correspond to the same topological 
order. 

To understand the above result, we notice 
that, from the structure of the fixed-point condi- 
tions, if (-F'fc;™^"^, F'f"''"^, is a solution, then 
i^pFl^'iJ^'"^ ,ripP^^'"^^ ,riAA^) is also a solution, where 
•qp = ±1, rjp — ±1, and rjA — ±1- However, the three 
phase factors rjp^ rjp, and tja do not lead to different 
fixed-point wave functions, since they only affect the 
total phase of the wave function and are unphysical. 

On the other hand, the fixed-point solutions may also 
contain phase factors that do correspond to different 
fixed-point wave functions. For example, the local or- 
thogonal transformation c"^^''^'" does not affect the fu- 
sion rule A^ijfc, where A/;^ is the number of links with 
|^o)-state and |/Q)-state. Such a local orthogonal trans- 
formation changes {Pkin^j'x ^ -^i'''"^ ^ generates a 
discrete family of the fixed-point wave functions. 

Similarly, we can consider the following local orthog- 
onal transformation \a) — ^ S/j ^i'^""''^'''' 1/^) ^^at acts 
on each vertex with states |«o), |jo) i l^o) on the three 
edges connecting to the vertex. Such a local orthogo- 
nal transformation also does not affect the fusion rule 
Nijk- The new local orthogonal transformation changes 

i^kin^'^x i^i'''"^ ^ ^'^d generates a family of the fixed- 
point wave functions parameterized by the orthogonal 
matrix 0^'T^''\ 

Now the question is that do those solutions related by 
local orthogonal transformations have the same time re- 
versal invariant topological order or not. We know that 
two gapped wave functions have the same time reversal 
invariant topological order if and only if they can be con- 
nected by local orthogonal transformation that connects 
to identity continuously. It is well known that an orthog- 
onal matrix whose determinant is — 1 does not connect to 
identity. Thus it appears that local orthogonal transfor- 
mations some times can generate different time reversal 
invariant topological orders. 

However, when we use the equivalent classes of local 
orthogonal transformations to define time reversal in- 
variant topological orders, we not only assume the local 
orthogonal transformations to connect to identity con- 
tinuously, we also assume that we can expand the lo- 
cal Hilbert spaces (say by increasing the range of the 
indices i and a that label the states on the edges and 
the vertices). The local orthogonal transformations can 
act on those enlarged Hilbert spaces and can connect to 
identity in those enlarged Hilbert spaces. Even when a 
local orthogonal transformation cannot be deformed into 
identity in the original Hilbert space, it can always be de- 
formed into identity continuously in an enlarged Hilbert 
space. Thus two real wave functions related by a local 
orthogonal transformation always have the same time re- 



versal invariant topological order. 

For example, an orthogonal matrix that acts 

on states |0) and |1) does not connect to identity within 
the space of two by two orthogonal matrices. However, 
we can embed the above orthogonal matrix into a three 
by three orthogonal matrix that acts on |0), |1), and |2): 
10 \ 

0—1 . Such a three by three orthogonal matrix 

-ly 

does connect to identity within the space of three by three 
orthogonal matrices. This completes our argument that 
all local orthogonal transformations can connect to iden- 
tity continuously at least in an enlarged Hilbert space. 

So, after factoring out the unimportant phase factors 
discussed above, the real solutions of the fixed-point con- 
ditions may uniquely correspond to time reversal invari- 
ant topological orders. The four types of real solutions 
discussed in the last section are examples of four different 
time reversal invariant topological orders. 

For the TV = 1 loop states and the = 1 string-net 
state, we only have two states on each link. In this case, 
we can treat the two states as the two states of an electron 
spin. The time reversal transformation now becomes 

t:co|0)+ci|l)^-ct|0)+c*|l) (70) 

on each link. For such a time reversal transformation 
f 2 = -1. The two iV 1 loop states and the TV = 1 
string-net state are not invariant under such an time- 
reversal transformation. 



XI. WAVE FUNCTION RENORMALIZATION 
FOR TENSOR PRODUCT STATES 

A. Motivation 

Once the fixed-point states have been identified and 
the labeling of topological orders has being found, we 
then face the next important issue: given a generic 
ground state wave function of a system, how to iden- 
tify the topological orders in the state? In other words, 

how to calculate the data {Nijk, Fkiri'-y'x'-^i'''"'^^ 
characterize the topological orders from a generic wave 
function? 

One way to address the above issue is to have a general 
renormalization procedure which flows other states in the 
same phase to the simple fixed-point state so that we can 
identify topological order from the resulting fixed-point 
state. That is, we want to find a local unitary transfor- 
mation which removes local entanglement and gets rid of 
unnecessary degrees of freedom from the state. How to 
find the appropriate unitary transformation for a specific 
state is then the central problem in this renormalization 
procedure. Such a procedure for one dimensional tensor 
product states(TPS) (also called matrix product states) 
has been given in Ref. 38. Here we will propose a method 



to renormalize two dimensional TPS, where nontrivial 
topological orders emerge. The basic idea is to use the 
gLU transformation discussed in section VII. Note that, 
through the gLU transformation, we can reduce the num- 
ber of labels in a region A to the minimal value without 
loosing any quantum information (see Fig. 6). This is 
because the gLU transformation is a lossless projection 
into the support space of the state in the region A. By 
performing such gLU transformations on overlapping re- 
gions repeatedly (see Fig. 2a), we can reduce a generic 
wave function to the simple fixed-point form discussed 
above. It should be noted that any state reducible in 
this way can be represented as a MERA'^^. 

In the following, we will present this renormalization 
procedure for two dimensional TPS where we find a 
method to calculate the proper gLU transformations. 
The tensor product states are many-body entangled 
quantum states described with local tensors. By mak- 
ing use of the entanglement information contained in the 
local tensors, we are able to come up with an efficient 
algorithm to renormalize two dimensional TPS. 

This algorithm can be very useful in the study of quan- 
tum phases. Due to the efficiency in representation, TPS 
has found wide application as variational ansatz states in 
the studies of quantum many-body systems^^"^^ . Sup- 
pose that in a variational study we have found a set of 
tensors which describe the ground state of a two dimen- 
sional many-body Hamiltonian and want to determine 
the phase this state belongs to. We can apply our renor- 
malization algorithm to this tensor product state, which 
removes local entanglement and flows the state to its 
fixed point. By identifying the kind of order present in 
the fixed-point state, we can obtain the phase informa- 
tion for the original state. 

In this section, we will give a detailed description of 
the algorithm and in the next section we will present its 
application to some simple (but nontrivial) cases. The 
states we are concerned with have translational symme- 
try and can be described with a translational invariant 
tensor network. To be specific, we discuss states on a 
hexagonal lattice. Generalization to other regular lat- 
tices is straight-forward. 
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FIG. 11: (Color online) Left: Tensor T representing a 2D 
quantum state on hexagonal lattice, i is the physical index, 
a, /3, 7 are inner indices. Right: a tensor product state where 
each vertex is associated with a tensor. The inner indices of 
the neighboring tensors connect according to the underlying 
hexagonal lattice. 

tensors by 

IV')- ti:T{T''T'\..T^...)\ni2...irn-) (71) 

where tTr denotes tensor contraction of all the connected 
inner indices on the links of the hexagonal lattice. 

A renormalization procedure of quantum states is com- 
posed of local unitary transformations and isometry maps 
such that the state fiows along the path IV''-^-'), 
|?/'^^)),...and finally towards a fixed point \'4>°°)- With the 
tensor product representation, flow of states corresponds 
to a flow of tensors T(°), T^^), T^^)... We wiU give the 
detailed procedure of how the tensors are mapped from 
one step to the next in the following section. 

C. Renormalization algorithm 

In one round of renormalization, we start from tensor 
T^") , do some operation to it which corresponds to local 
unitary transformations on the state, and map T^"^ to 
y(n+i)^ The whole procedure can be broken into two 
parts: the F-move and the P move, in accordance with 
the two steps introduced in the previous section. 



B. Tensor product states 

Consider a two-dimensional spin model on a hexagonal 
lattice with one spin (or one qudit) living at each vertex. 
The Hilbert space of each spin is D-dimensional. The 
state can be represented by assigning to every vertex a 
set of tensors T^^^, where i labels the local physical di- 
mension and takes value from 1 to D. a,/3,7 are inner 
indices along the three directions in the hexagonal lat- 
tice respectively. The dimension of the inner indices is d. 
Note that the figures in this note are all sideviews with 
inner indices in the horizontal plane and the physical in- 
dices pointing in the vertical direction, if not specified 
otherwise. The wave function is given in terms of these 



1. Step 1: F-move 

In the F-move, we take a ^—^^ configuration in the ten- 
sor network and map it to a configuration by doing a 

local unitary operation. We will see that the tensor prod- 
uct representation of a state leads to a natural way of 
choosing an appropriate unitary operation for the renor- 
malization of the state. 

In order to do so, first we define the double tensor T 
of tensor T as 

Ta'/3'7',Q/J7 = '^{Ta'p'y)* X T*^^ (72) 
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FIG. 12: (Color online) Double tensor T represented as two 
layers of tensor T with the physical indices contracted. The 
gray layer is the lower layer. 




FIG. 13: (Color online) F-move in the renormalization proce- 
dure: (1) combining double tensors Ti and T2 on neighboring 
sites into a single double tensor T (2)Splitting double tensor T 
into tensor T (3) SVD decomposition of tensor T into tensors 
Ta and Tb. 



Graphically the double tensor T is represented by two 
layers of tensor T with the physical indices connected. 
The tensor T giving rise to the same double tensor T is 
not unique. Any tensor T' which differs from T by an 
unitary transformation U on physical index i gives the 
same T as U and cancels out in the contraction of i. 
On the other hand, an unitary transformation on i is the 
only degree of freedom possible, i.e. any T' which gives 
rise to the same T as T differs from T by a unitary on 
i. Therefore, in the process of turning a tensor T into a 
double tensor T and then split it again into a different 
tensor T', we apply a non-trivial local unitary operation 
on the corresponding state. A well designed way of split- 
ting the double tensor will give us the appropriate unitary 
transformation we need, as we show below. 

F-move has the following steps. First, construct dou- 
ble tensors for two neighboring sites on the lattice and 
combine them into a single double tensor with 8 inner 
indices. 

a' P' 5' t' ,aP5<i = '^^'^ l,a' fj' -j' ,a[j-t X ^^2,5' e'f' .St-y (73) 

7', 7 

Note that with respect to the bipartition of indices 



a' jj'S'e' and a/3Se, T is Hermitian 

'^a'l3'S'e',al3Se — aPSe,a' jS' S' t')* (74) 

and positive semidefinite. Therefore it has a spectral 
decomposition with positive eigenvalues {A^ > 0}. The 
corresponding eigenvectors are {T-^} 

j 

This spectral decomposition lead to a special way of de- 
composing double tensor T into tensors. Define a rank 8 
tensor r(as shown in Fig. 13 after step 2) as follows: 

Ti"?: = E Va; {fLnr) * X Ti,s. (76) 

i 

T has four inner indices a, /?, 5, e of dimension d and four 
physical indices l,m,n,r also of dimension d which are 
in the direction of a,(3,S,e respectively. As {T-*} form 
an orthonormal set, it is easy to check that T gives rise 
to double tensor T. Going from Ti and T2 to T, we have 
implemented a local unitary transformation on the phys- 
ical degrees of freedom on the two sites, so that in T the 
physical indices and the inner indices represent the same 
configuration. In some sense, we are keeping only the 
physical degrees of freedom necessary for entanglement 
with the rest of the system while getting rid of those 
that are only entangled within this local region. Now we 
do a singular value decomposition of tensor T in the di- 
rection orthogonal to the link between Ti and T2 and T 
is decomposed into tensors Ta and Tf,. 

rplmnr \ ^ rpln ^ rpmr {'7'7\ 

X 

This step completes the F-move. Ideally, this step should 
be done exactly so we are only applying local unitary 
operations to the state. Numerically, we keep some large 
but finite cutoff dimension for the SVD step, so this step 
is approximate. 

On a hexagonal lattice, we do F-move on the chosen 
neighboring pairs of sites (dash-circled in Fig. 14), so 
that the tensor network is changed into a configuration 
shown by thick dark lines in Fig. 14. Physical indices 
are omitted from this figure. Now by grouping together 
the three tensors that meet at a triangle, we can map 
the tensor network back into a hexagonal lattice, with 
1/3 the number of sites in the original lattice. This is 
achieved by the P-move introduced in the next section. 

2. Step 2: P-move 

Now we contract the three tensors that meet at a tri- 
angle together to form a new tensor in the renormalized 
lattice as shown in the first step in Fig. 15 

rpl _ rpnr rplm rppq fjo\ 
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FIG. 14: Original hexagonal lattice (gray line) and renormal- 
ized lattice (black line) after F-move has been applied to the 
neighboring pairs of sites circled by dash line. 



Ta 




FIG. 15: (Color online) P-move in the renormalization pro- 
cedure: (1) contracting three tensors that meet at a triangle 
Ta,Ti,Tc to form a new tensor T'^' on one site of the renor- 
malized hexagonal lattice. (2)constructing the double tensor 
T'^^ from T^^' so that we can start to do F-moves again. 

where / is the physical index of the new tensor which in- 
cludes all the physical indices of T^, T^, T^: I, m, n, r, p, q. 
Note that in the contraction, only inner indices are con- 
tracted and the physical indices are simply group to- 
gether. 

Constructing the double tensor T from T, we get the 
renormalized double tensor on the new hexagonal lattice 
which is in the same form as Ti,T2 and we can go back 
again and do the F-move. 



3. Complications: corner double line 

One problem with the above renormalization algorithm 
is that, instead of having one isolated fixed-point tensor 
for each phase, the algorithm has a continuous family 
of fixed points which all correspond to the same phase. 
Consider a tensor with structure shown in Fig. 16. The 
tensor is a tensor product of three parts which include 
indices {a, /3, ii}, {7, iS, ^2}, {e, A, 13} respectively. It can 
be shown that this structure remains invariant under the 
renormalization flow. Therefore, any tensor of this struc- 
ture is a fixed point of our renormalization flow. However, 
it is easy to see that the state it represents is a tensor 
product of loops around each plaquette, which can be 
disentangled locally into a trivial product state. There- 
fore, the states all have only short-range entanglement 
and correspond to the topologically trivial phase. The 



FIG. 16: (Color online) The corner double line tensor which 
is a fixed point of the renormalization algorithm. The three 
groups of indices {a, /3, ii}, {7, (5, 12}, {e, A, is}, are entangled 
within each group but not between the groups. 

trivial phase has then a continuous family of fixed-point 
tensors. This situation is very similar to that discussed in 
Ref. 37,69. We will keep the terminology and call such a 
tensor a corner double line tensor. Not only does corner 
double line tensor complicate the situation in the trivial 
phase, it leads to a continuous family of fixed points in 
every phase. It can be checked that the tensor product 
of a corner double line with any other fixed-point tensor 
is still a fixed-point tensor. The states they correspond 
to differ only by small loops around each plaquette and 
represent the same topological order. Therefore any sin- 
gle fixed-point tensor gets complicated into a continuous 
class of fixed-point tensors. In practical application of the 
renormalization algorithm, in order to identify the topo- 
logical order of the fixed-point tensor, we need to get rid 
of such corner double line structures. Due to their sim- 
ple structure, this can always be done, as discussed in 
the next section. 



XII. APPLICATIONS OF THE 
RENORMALIZATION FOR TENSOR PRODUCT 
STATES 

Now we present some examples where our algorithm 
is used to determine the phase of a tensor product state. 
The algorithm can be applied both to symmetry breaking 
phases and topological ordered phases. In the study of 
symmetry breaking/topological ordered phases, suppose 
that we have obtained some tensor product description 
of the ground state of the system Hamiltonian. We can 
then apply our algorithm to the tensors, flow them to the 
fixed point, and see whether they represent a state in the 
symmetry breaking phase/topological ordered phase or a 
trivial phase. 

For system with symmetry/topological order related to 
gauge symmetry, it is very important to keep the sym- 
metry/gauge symmetry in the variational approach to 
ground state and search within the set of tensors that 
have this symmetry /gauge symmetry^''. The resulting 
tensor will be invariant under such symmetries/gauge 
symmetries, but the state they correspond to may have 
different orders. In the symmetry breaking case, the 
state could have this symmetry or could spontaneous 
break it. In the topological ordered phase, the state 
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FIG. 17: (Color online) Renormalization procedure on square 
lattice part 1: (1) decomposing double tensor T into tensor 
T (2) SVD decomposition of T in two different directions, 
resulting in tensors Ta, Th and Tc, Td respectively. 



could have nontrivial topological order or be just triv- 
ial. Our algorithm can then be applied to decide which 
is the case. In order to correctly determine the phase for 
such symmetric tensors, it is crucial that we maintain the 
symmetry/gauge symmetry of the tensor throughout our 
renormalization process. We will discuss in detail two 
cases: the Ising symmetry breaking phase and the Z2 
topological ordered phase. For simplicity of discussion 
and to demonstrate the generality of our renormalization 
scheme, we will first introduce the square lattice version 
of the algorithm. All subsequent applications are carried 
out on square lattice. (Algorithm on a hexagonal lattice 
would give qualitatively similar result, though quantita- 
tively they might differ, e.g. on the position of critical 
point.) 



A. Renormalization on square lattice 

Tensor product states on a square lattice are repre- 
sented with one tensor T^^^^g on each vertex, where i is 
the physical index and a 0^6 are the four inner indices 
in the up, down, left, right direction respectively. We will 
assume translational invariance and require the tensor to 
be the same on every vertex. The renormalization pro- 
cedure is be a local unitary transformation on the state 
which flows the form of the tensor until it reaches the 
fixed point. It is implemented in the following steps. 
First, we form the double tensor T from tensor T 

i 

Then do the spectral decomposition of positive operator 
T into 

j 




FIG. 18: Original square lattice (gray line) and renormalized 
lattice (black line) after the operations in Fig. 17 has been 
applied. Physical indices of the tensors are not drawn here. 




FIG. 19: (Color online) Renormalization procedure on square 
lattice part 2: (1) combining four tensors that meet at a 
square into a single tensor (2) constructing a double tensor 
from it. 

and form a new tensor 3^™"^ 

j 

l,m,n,r are physical indices in the up, down, left, right 
directions respectively. This is illustrated in step 1 of 
Fig. 17. This step is very similar to the second step in 
the F-move on hexagonal lattice. Next we do SVD de- 
composition of tensor T. For vertices in sublattice A we 
decompose between the up-right and down-left direction 
as shown in step 2 of Fig. 17. For vertices in sublat- 
tice B we decompose between the up-left and down-right 
direction as shown in step 3 of Fig. 17. 

rplmnr \ ^ rplr ^ rpmn 

A 

rplmnr \ ^ rplr ^ rpmn 

^a^jS — ^^c,aSX ^ -'-b,l3f\ 

A 

After the decomposition, the original lattice (gray lines 
in Fig. 18) is transformed into the configuration shown 
by thick dark lines in Fig. 18. Physical indices are omit- 
ted from this figure. If we now shrink the small squares, 
we get a tensor product state on a renormalized square 
lattice. Fig. 19 shows how this is done. 

In step 1, we contract the four tensors that meet at a 
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small square 

rpl _ \ " rpmirii rpm-: 
^ap-ti ^ a,l}ux b,a 



where / stands for the combination of all physical indices 
niijirii, i = 1...4. In step 2, we construct a double tensor 
T form tensor T and completes one round of renormal- 
ization. Now we can go back to step 1 in Fig. 17 and 
flow the tensor further. 

Now we are ready to discuss two particular examples, 
the Ising symmetry breaking phase and the Z2 topologi- 
cal ordered phase, to demonstrate how our algorithm can 
be used to determine the phase of a tensor product state. 



B. Ising symmetry breaking phase 

A typical example for symmetry breaking phase transi- 
tion is the transverse field Ising model. Consider a square 
lattice with one spin 1/2 on each site. The transverse field 
Ising model is 



H 



I sing 



where {ij} are nearest neighbor sites. The hamiltonian 
is invariant under spin flip transformation J^j, Xk for any 
e. 

When e = 0, the ground state spontaneously breaks 
this symmetry into either the all spin up state |00...0) or 
the all spin down state |ll...l). In this case any global 
superposition a|00...0) -|-/3|ll...l) represents a degenerate 
ground state. When e = 00, the ground state has all spin 
polarized in the X direction (| + +...+)) and does not 
break this symmetry. 

In the variational study of this system, we can require 
that the variational ground state always have this sym- 
metry, regardless if the system is in the symmetry break- 
ing phase or not. Then we will find for e = the ground 
state to be |00...0) -I- |ll...l). Such a global superposi- 
tion represents the spontaneous symmetry breaking. For 
e = cx), we will find the ground state to be | -I- -f ...-!-) and 
does not break the symmetry. For < e < cx), we will 
need to decide which of the previous two cases it belongs 
to. We can first find a tensor product representation of 
an approximate ground state which is symmetric under 
the spin flip transformation, then apply the renormaliza- 
tion algorithm to find the fixed point and decide which 
phase the state belongs to. Below we will assume a sim- 
ple form of tensor and demonstrate how the algorithm 
works. 

Suppose that the tensors obtained from the variational 
study TJ^p-ygj where i,a,f3,j,S can be or 1, takes the 
following form 



7^0 _ \Q+/3+7+<5 

rpl ^ ,4_(Q + /3+^ + i) 




FIG. 20: Quantity X2/ X\ obtained by taking the ratio of the 
contraction value of the double tensor in two different ways. 

is invariant under change of scale, basis transforma- 
tion and corner double line structures of the double tensor 
and can be used to distinguish different fixed point tensors. 
For clarity, only one layer of the double tensor is shown. The 
other layer connects in exactly the same way. 



A is a parameter between and 1. Under an X operation 
to the physical index, the tensor is changed to T 



rfl 



(79) 



T can be mapped back to T by switching the 0,1 label 
for the four inner indices aPjS. Such a change of basis 
for the inner indices does not change the contraction re- 
sult of the tensor and hence the state that is represented. 
Therefore, the state is invariant under the spin flip trans- 
formation Y[k Xk and we will say that the tensor has this 
symmetry also. 

When A = 0, the tensor represents state |00...0) -I- 
|ll...l), which corresponds to the spontaneous symme- 
try breaking phase. We note that the A = tensor is 
a direct sum of dimension- 1 tensors. Such a direct-sum 
structure corresponds to spontaneous symmetry break- 
ing, as discussed in detail in Ref. 37. When A = 1, the 
tensor represents state | + +...+) which corresponds to 
the symmetric phase. When < A < 1, there must be 
a phase transition between the two phases. However, as 
A goes from to 1, the tensor varies smoothly with well 
defined symmetry. It is hard to identify the phase tran- 
sition point. Now we can apply our algorithm to the 
tensor. First, we notice that at A = or 1, the tensor is 
a fixed point for our algorithm. Next, we find that for 
A < 0.358, the tensor flows to the form with A = 0, while 
for A > 0.359, it flows to the form with A = 1. Therefore, 
we can clearly identify the phase a state belongs to using 
this algorithm and find the phase transition point. 

Note that in our algorithm, we explicitly keep the spin 
flip symmetry in the tensor. That is, after each renormal- 
ization step, we make sure that the renormalized tensor 
is invariant under spin flip operations up to change of 
basis for the inner indices. If the symmetry is not care- 
fully preserved, we will not be able to tell the two phases 
apart. 

We also need to mention that for arbitrary A, the fixed 
point that the tensor flows to can be different from the 
tensor at A = or 1 by a corner double line structure. We 
need to get rid of the corner double line structure in the 
result to identify the real fixed point. This is possible 
by carefully examining the fixed point structure. An- 
other way to distinguish the different fixed points with- 
out worrying about corner double lines is to calculate 
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FIG. 21: (Color online) X2/X1 for tensors (79) under the 
renormalization flow. As the number of RG steps increases, 
the transition in X2/X1 becomes sharper and finally ap- 
proaches a step function at fixed point. The critical point 
is at Ac ~ 0.358. 

some quantities from the fixed-point tensors that are in- 
variant with the addition of corner double lines. We 
also want the quantity to be invariant under some triv- 
ial changes to the fixed point, such as a change in scale 
T — rjT or the change of basis for physical and inner 
indices. One such quantity is given by the ratio of X2 

and Xi defined as Xi = {j2a'Y .a-y'^a'a'Yj' -aa-ry^ , and 

^2 = J2a'0'YS',al3y5^a'a'YS',aajS X T /^^ p> S' j' ,13 PS^ ■ Fig. 

20 gives a graphical representation of these two quanti- 
ties. In this figure, only one layer of the double tensor is 
shown. The other layer connects in the exactly the same 
way. It is easy to verify that X2 / Xi is invariant under the 
change of scale, basis transformation and corner double 
lines. 

We calculate X2/X1 along the renormalization flow. 
The result is shown in Fig. 21. At the A = fixed 
point, X2IX1 = 0.5 while at A = 1, X2/X1 = 1. As we 
increase the number of renormalization steps, the transi- 
tion between the two fixed points becomes sharper and 
finally approaches a step function with critical point at 
Ac ~ 0.358. Tensors with A < Ac belongs to the symme- 
try breaking phase while tensors with A > Ac belongs to 
the symmetric phase. 

C. Z2 topological ordered phase 

The algorithm can also be used to study topological 
order of quantum states. In this section, we will demon- 
strate how the algorithm works with Z2 topological order. 

Consider again a square lattice but now with one spin 
1/2 per each link. A simple Hamiltonian on this lattice 
with Z2 topological order can be defined as 



where p means plaquettes and i G p is all the spin l/2s 
around the plaquette and v means vertices and j e u is 
all the spin l/2s connected to the vertex. The ground 
state wave function of this Hamiltonian is a fixed-point 
wave function and corresponds to the — 1 loop state 
with 77 = 1 as discussed in the previous section. 

The ground state wave function has a simple tensor 
product representation. For simplicity of discussion we 
split every spin 1/2 into two and associate every vertex 
with four spins. The tensor T^^^g has four physical 
indices i, j, kj = 0,1 and three inner indices a, l3,j,S = 
0,1. 

T:^mm^^, ifmod(i+j-ffc-H;,2) = 
all other terms being 

It can be checked that T^^ is a fixed-point tensor of our 
algorithm. This tensor has a Z2 gauge symmetry. If we 
apply Z operation to all the inner indices, where Z maps 
to and 1 to —1, the tensor remains invariant as only 
even configurations of the inner indices are nonzero in 
the tensor. 

Consider then the following set of tensor parameterized 
by g 

Tl^H = 5'+^+'^+', if mod(z +j+k + l,2) = Q 
all other terms being 

Ax g — 1, this is exactly T^^ and the corresponding state 
has topological order. At 5 = 0, the tensor represents 
a product state of all and we denote the tensor as Tq. 
At some critical point in g, the state must go through 
a phase transition. This set of tensors are all invariant 
under gauge transformation ZZZZ on their inner indices 
and the tensor seems to vary smoothly with g. One way 
to detect the phase transition is to apply our algorithm. 
We find that, at 5 > ^c, the tensors flow to Tg^, while 
at g < gc, the tensors flow to Tq. We determine gc to be 
between 0.804 ~ 0.805. As this model is mathematically 
equivalent to two dimensional classical Ising model where 
the transition point is known to great accuracy, we com- 
pare our result to that result and flnd our result to be 
within 1% accuracy {gc = 0.8022). Again in the renor- 
malization algorithm, we need to carefully preserve the 
Z2 gauge symmetry of the tensor so that we can correctly 
determine the phase of the states. 

The fixed-point tensor structure might also be compli- 
cated by corner double line structures, but it is always 
possible to identify and get rid of them. Similarly, we 
can calculate the invariance quantity X2/X1 to distin- 
guish the two fixed points. X2/X1 — 1 for while 
X2/X1 = 0.5 for To. The resufl is plotted in Fig. 22 and 
we can see that the transition in X2/X1 approaches a 
step function after a large number of steps of RG, i.e. at 
the fixed point. The critical point is at gc = 0.804. For 
g < gc, the tensor belongs to the trivial phase, while for 
g > gc, the tensor belongs to the Z2 topological ordered 
phase. 
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FIG. 22: (Color online) X2/X1 for tensors (80) under the 
renormalization flow. As the number of RG steps increases, 
the transition in X2/X1 becomes sharper and finally ap- 
proaches a step function at fixed point. The critical point 
is at Qc = 0.804. 



Our algorithm can also be used to demonstrate the 
stability of topological order against local perturbation. 
As is shown in Ref. 70, local perturbations to the Z2 
Hamiltonian correspond to variations in tensor that do 
not break the Z2 gauge symmetry. We picked tensors 
in the neighborhood of T^^ which preserve this gauge 
symmetry randomly and applied our renormalization al- 
gorithm (gauge symmetry is kept throughout the renor- 
malization process). We find that as long as the variation 
is small enough, the tensor fiows back to Ti^ , up to a cor- 
ner double line structure. This result shows that the Z2 
topological ordered phase is stable against local pertur- 
bations. 



XIII. SUMMARY 

In this paper, we discuss a defining relation between 
local unitary transformation and quantum phases. We 
argue that two gapped states are related by a local uni- 
tary transformation if and only if the two states belong 
to the same quantum phase. 

We can use the equivalent classes of local unitary trans- 
formations to define "patterns of long range entangle- 
ment" . So the patterns of long range entanglement cor- 
respond to universality classes of quantum phases, and 
are the essence of topological orders. '^^ 

As an application of this point of view of quantum 
phases and topological order, we use the generalized lo- 
cal unitary transformations to generate a wave function 
renormalization, where the wave functions flow within 
the same universality class of a quantum phase (or the 
same equivalent class of the local unitary transforma- 
tions). In other words, the renormalization flow of a wave 
function does not change its topological order. 

Such a wave function renormalization allows us to clas- 



sify topological orders, by classifying the fixed-point wave 
functions and the associated fixed-point local unitary 
transformations. 

First, we find that the fixed-point local uni- 
tary transformations are described by the data 
(%^.^^feLT/'^'''"'') that satisfy 

N 



N 



N 



J2 N3^rn'Nk„^l' ^J^^kjn'Ni 



l*ni ) 



(81) 
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^mcx^^m' a' p' 



^kln,xS \-^kln,xS > 

Nkjt' Nu-n' Nit 

EST^ pii'm,af} -pitn.ipxpjkt, 

/ J / J / J knt.T]!^ Ips.K^ lsq.6 
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\0F \ ^ Tpm,kn,Px-p' 



qps,<p'-y ' 



(82) 



jkj.aP 



Nkjk* 

^klk.xS fc* 



for all fc, z, / satisfying Nku- > 0. 



(83) 



[see eqns. (9, 30, 47 12, 18, 23, 33, 36)]. From the data 
{Nijk, Fkin '^s^ ,Pi'''°'^) we can further find out the fixed- 
point wave function by solving the following equations 



for A\i = 0,...,iV: 



A' 



m, A,7 

ikj^aj3 kji,aP^ 

'i>lj,c.p^O, ifA^.fe, =0, 
deti<^'kj,o.f}) + 0. 



(84) 



[see eqns. (40, 38, 42, 44, 45, 46)]. 



The combined data (A^^jfc, Ffc^"^"/, -P/"''"^, ^4*) that 
satisfy the conditions (81), (82), (83) and (84) classify 
a large class of topological orders. We see that the prob- 
lem of classifying a large class of topological orders be- 
comes the problem of solving a set of non-linear algebraic 
equations (81), (82), (83), and (84). 
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The combined data {N,jk,F;^i:''^f ^P^^'^^A^) that 
satisfy the conditions (81), (82), (83), and (84) also clas- 
sify a large class of time reversal invariant topological 
orders, if we restrict ourselves to real solutions. The so- 
lutions related by local orthogonal transformations all be- 
long to the same phase, since the local orthogonal trans- 
formations always connect to identity if we enlarge the 
Hilbert space. 

We like to point out that we cannot claim that the 
solutions (%fc,i^fer;,f ,i^'^'"'',A») classify ah topolog- 
ical orders since we have assumed that the fixed-point 
local unitary transformations are described by tensors 
of finite dimensions. It appears that chiral topological 
orders, such as quantum Hall states, are described by 
tensors of infinite dimensions. 

We note that the data (%fe, i^^f™^"/ , i^'^^'^"^, just 
characterize different fixed-point wave functions. It is 
not guaranteed that the different data will represent dif- 
ferent topological orders. However, for the simple so- 
lutions discussed in this paper, they all coincide with 
string-net states, where the topological properties, such 
as the ground state degeneracy, number of quasiparti- 
cle types, the quasiparticle statistics etc, were calculated 
from the data. From those topological properties, we 
know that those different simple solutions represent dif- 
ferent topological orders. Also (^fc;™^^^, -Pf"''"'^, ^*) are 
real for the simple solutions discussed here. Thus they 
also represent topological orders with time reversal sym- 
metry. (Certainly, at the same time, they represent stable 
topological orders even without time reversal symmetry.) 

We also like to point out that our description of fixed- 
point local unitary transformations is very similar to the 
description of string-net states. However, the conditions 
(81), (82), (83), and (84) on the data appear to be weaker 
than (or equivalent to) those^^ on the string-net data. 
So the fixed-point wave functions discussed in this paper 
may include all the string-net states (in 2D). 

Last, we present a wave function renormalization 
scheme, based on the gLU transformations for generic 
TPS. Such a wave function renormalization always flows 
within the same phase (or within the same equivalence 
class of LU transformations). It allows us to determine 
which phase a generic TPS belongs to by studying the re- 
sulting fixed-point wave functions. We demonstrated the 
effectiveness of our method for both symmetry break- 
ing phases and topological ordered phases. We find that 
we can even use tensors that do not break symmetry 
to describe spontaneous symmetry breaking states: if a 
state described by a symmetric tensor has a spontaneous 
symmetry breaking, the symmetric tensor will fiow to a 
fixed-point tensor that has a form of direct sum. 
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Appendix: Equivalence relation between quantum 
states in the same phase 

In section III and section IV, we argued about the 
equivalence relation between gapped quantum ground 
states in the same phase. We concluded that two states 
are in the same phase if and only if they can be connected 
by local unitary evolution or constant depth quantum cir- 
cuit: 

|$(1)) - |$(0)) iff |$(1)) =r[e-*-^" '^5^(3)]|$(0)) 

i$(i)) ^ i$(o)) iff mi)) = u^jm) 

Now we want to make these arguments more rigorous, 
by stating clearly what is proved and what is conjectured, 
and by giving precise definition of two states being the 
same, the locality of operators, etc. We will show the 
equivalence in the following steps: 

1. if two gapped ground states are in the same phase, 
then they are connected by local unitary evolution 

2. a gapped ground state remains in the same phase 
under local unitary evolution 

3. a local unitary evolution can be simulated by a 
constant depth quantum circuit and vice-verse. 

(All these discussions can be generalized to the case 
where the system has certain symmetries. H{g) and Udrc 
used in the equivalence relation will then have the same 
symmetry as the system Hamiltonian H.) 

First, according to the definition in section III, two 
states |$(0)) and |$(1)) are in the same phase if we can 
find a family of local Hamiltonians H{g), g e [0, 1] with 
I $(5)) being its ground state such that the ground state 
average of any local operator O, {^{g)\0\^{g)) changes 
smoothly from g = to g = 1. Here we allow a more gen- 
eral notion of locality for the Hamiltonian'^^ and require 
H{g) to be a sum of local operators Hz{g)- 

H{g) = ^^(5) (85) 

Z<£Z 

where Hz{g) is a Hermitian operator defined on a com- 
pact region Z. X]ze2 sums over a set Z of regions. The 
set Z contains regions that differ by translations. The 
set Z also contains regions with different sizes. However, 
Hz{g) approaches zero exponentially as the size of the 
region Z approaches infinity. Or more precisely, for all 
sites u in the lattice 

\\Hz{g)\\ \Z\ expi^idiamiZ)) - 0(1) (86) 

zez,Z3u 

for some positive constant /i. Here X^zez Z3u sums over 
all regions in the set Z that cover the site u, \\...\\ denotes 
operator norm, \Z\ is the cardinality of Z, and diam(Z) 
is the diameter of Z. Therefore, instead of being exactly 
zero outside of a finite region, the interaction terms can 
have an exponentially decaying tail. 

If |<i>(0)) and |$(1)) are gapped ground states of H{0) 
and H{1), then we assume that for all {^{g)\0\^{g)) to 
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be smooth, H(g) must remain gapped for all g. If H{g) 
closes gap for some gc, then there must exist a local op- 
erator O such that {^{g)\0\^{g)) has a singularity at 
gc- We call the gapped H{g) an adiabatic connection 
between two states in the same phase. 

The existence of an adiabatic connection gives rise to 
a local unitary evolution between |<&(0)) and [^(l)). A 
slight modification of Lemma 7.1 in Ref. 33 gives that 

Theorem 1 Let H{g) be a differentiable family of lo- 
cal Hamiltonians and \^{g)) be its ground state. If 
the excitation gap above \^{g)) is larger than some fi- 
nite value A for all g, then we can define H{g) — 
i J dtF{t)exp{iH{g)t) {dgH{g)) exp{—iH{g)t), such that 

|$(1)> =r[e-*/o 'i9-f^(5)]|$(0)). 

where F{t) is a function which has the following proper- 
ties: 1. The fourier transform of F{t), F{uj) is equal to 
— l/ui for |a;| > A. 2. F{uj) is infinitely differentiable. 
3. F{t) = —F{~t). Under this definition, H{g) is local 
(almost) and satisfies 



Hz{9).0b </i'(dist(Z,B))|Z| Hz{g) 



\Of 



(87) 

where Ob is any operator supported on site B, dist(Z, B) 
is the smallest distance between B and any site in Z , and 
h' [I) is a function which decays faster than any negative 
power of I. [...] denotes commutator of two operators. 
This is called the quasi-adiabatic continuation of states. 

Therefore, we can show that if |$(0)) and |$(1)) are 
in the same phase, then we can find a local (as de- 
fined in (87)) Hamiltonian Hz{g) such that |<i>(l)) = 
7-[e"*/o '^9'f^(f)]|$(0)). In other words, states in the same 
phase are equivalent under local unitary evolution. 

Note that here we map |$(0)) exactly to |<I'(1)). There 
is another version of quasi-adiabatic continuation''^, 
where the mapping is approximate. In that case, Hz{g) 
can be defined to have only exponentially small tail out- 
side of a finite region instead of a tail which decays faster 
than any negative power. T[e~^^o ''^■^(^^]|$(0)) will not 
be exactly the same as [^(l)), but any local measurement 
on them will give approximately the same result. 

Next we want to show that the reverse is also true. 
Suppose that |<&(0)) is the gapped ground state of a 
local Hamiltonian -ff(O), -ff(O) = J2z ^z{0) and each 
Hz{0) is supported on a finite region Z. Apply a lo- 
cal unitary evolution generated by H{s) to |$(0)) and 
take it to |$(g)), 5 = - 1, |$(g)) = UgMO)), where 



where h{l) decays faster than any negative power of 
Therefore, Hz{g) remains local up to a tail which decays 
faster than any negative power. H{g) then forms a local 
gapped adiabatic connection between |'&(0)) and | $(!)). 

To detect for phase transition, we must check whether 
the ground state average value of any local operator O, 
(<i>(g)|0|$(.g)), has a singularity or not. (<i>(g)|0|$(.g)) 
($(0)|;7g^OC/g|$(0)). Using the Lieb-Robinson bound 
given in Ref. 33, we find U^OUg remains local (in the 
sense of (87)) and evolves smoothly with g. Therefore, 
{^{g)\0\^{g)) changes smoothly. In fact, because H{g) is 
differentiable, the derivative of (<I>(g)|0|$(g)) to any or- 
der always exists. Therefore, there is no singularity in the 
ground state average value of any local operator O and 
|$(0)) and |<I'(1)) are in the same phase. We have hence 
shown that states connected with local unitary evolution 
are in the same phase. 

This completes the equivalence relation stated in (1). 
The definition of locality is slightly different in different 
cases, so there is still some gap in the equivalence rela- 
tion. However, we believe that the equivalence relation 
should be valid with a slight generalization of the defini- 
tion of locality. 

Lastly, we want to show that the equivalence relation 
is still valid if we use constant depth quantum circuit 
instead of local unitary evolution. This is true because 
we can always simulate a local unitary evolution using a 
constant depth quantum circuit and vice-verse. 

To simulate a local unitary evolution of the form 
T[e^^-^o '^3 H(g)-^^ gj.g|. (divide the total time into small seg- 

set of local opera- 



(m+i)at 



ments 5t such that T[e '^^^ 
H{mSt) = Y^z^zimSt). The 
tors {Hz{m5t)} can always be divided into a finite 
number of subsets {/J^j (mJt)}, {iJ^i (m^t)}... such 
that elements in the same subset commute with each 
other. Then we can do Trotter expansion and approx- 



imate e 



-iStH{mSt) 



as e 

-iStH j n 



-iStH(mSt) 



(C/(i)f/(2)...)", where 
L/d) = n^e'"'"-l^", f/(2) = n.e"''*''-^/"... Each 
term in J/^^ C/*^^-* commute, therefore they can be im- 
plement as a piece-wise local unitary operator. Putting 
these piece-wise local unitary operators together, we have 
a quantum circuit which simulates the local unitary evo- 
lution. The depth of circuit is proportional to n x {l/5t). 
It can be shown that to achieve a simulation with con- 
stant error, a constant depth circuit would suffice^^. 
On the other hand, to simulate a constant depth 



quantum circuit 

(fe) 



pwl pwl ' ' 



pwl 



where U, 



(fe) 



pwl 



time dependent Hamiltonian as H{t) = X^i-^i'^^ such 



that for {k - l)5t < t < kSt, e'^'^"'^^ = [/f It is easy 
to check that r[e"'-/'o""'*f ^f(ff)] = jj^irc- The simulation 
time needed is M5t. We can always choose a finite 5t such 



Ug = T[e-'iS dsH(s)Y is then ground state of 11. Ul with a focal unitary evolution, we can define the 

H{g) = UgH{0)Ul = TzUaHziml = TzHzig)- 
Under unitary transformation Ug the spectrum of H(Q) 
doesn't change, therefore H{g) remains gapped. To show 
that H{g) also remains local, we use the Lieb-Robinson 
bound derived in Ref. 33, which gives 



that 



\[Hz{g).OB]\\ - \\[UgHz{mlOB]\\ 

< h{dist{Z,B))\Z\\\Hzm\ 



H, 



(fc) 



is finite and the evolution time is finite. 



\Oi 



We have shown that constant depth quantum circuit 
and local unitary evolution can simulate each other. The 
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equivalence relation (1) can therefore also be stated in terms of constant depth circuit as in (4). 
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The name "topological order" is motivated by the low en- 
ergy effective theory of the chiral spin states, which is a 
topological quantum field theory/'*. 

The conjecture is for symmetry breaking states with finite 
energy gaps. However, it is unclear whether it is possible to 
find a definition of "short range quantum entanglement" to 
make the conjecture valid for gapless symmetry breaking 
states. 

For each Hamiltonian, we can introduce an algebra of local 



operators. If the energy gap for H{g) is finite for all g in the 
segment [0, 1], Ref. 34 show that the algebra of local oper- 
ators for different g's are isomorphic to each other at low 
energies, using a quasi-adiabatic continuation of quantum 
states. 

A state is called a gapped state if there exist a local Hamil- 
tonian such that the state is the gapped ground state of 
the Hamiltonian. 



